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SUMMARY 


This  report  discusses  the  data  needed  to  perform  radiative  transfer 
calculations  in  nonhydrogenic  gases  in  local  thermodynamic  equilibrium  and 
presents  some  approximate  methods  for  computing  the  radiative  energy 
transferred  by  spectral  lines  where  the  properties  of  the  gas  are  uniform. 

The  methods  currently  available  for  calculating  the  cross  sections  of 
radiative  processes  are  described  and  compared.  An  accurate  method  for 
calculating  the  species  composition  of  nitrogen  is  described  and  the  results 
of  such  a  calculation  are  presented.  The  important  line  broadening 
mechanisms  are  discussed  and  the  potentially  accurate,  modern  theories  of 
line  broadening  are  outlined.  The  results  of  these  theories  are  used  to 
justify  approximate  line  profiles  which  are  simple  enough  for  use  in 
radiative  transfer  calculations. 

Simple  approximations  to  the  exact  curves  of  growth  of  intensity  are 
described  for  lines  with  Doppler  profiles  and  for  lines  with  profiles  of 
a  class  which  includes  the  dispersion  and  quasi-static  forms.  The  concept 
of  the  effective  width  of  a  line  intensity  profile  is  introduced  and 
techniques  are  developed  for  dealing  with  the  overlapping  of  the  intensity 
profiles  of  small  groups  of  closely  spaced  lines  (as,  for  example,  in  a 
multiplet ) . 
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general  introduction  and  outline  or  REPORTS 

The  transfer  of  energy  by  radiative  processes  in  gases  has  been  a 
subject  of  study  for  many  years  by  astrophysicists  and  spectroscopists. 

The  aim  of  these  workers  has  for  the  most  part  been  to  examine  the  spec¬ 
tral  distribution  of  radiation  and  to  use  it  to  predict  the  internal 
structure  and  state  of  the  gas  from  which  the  radiation  emerged.  More 
recently,  the  high  temperatures  associated  with  nuclear  explosions  and 
with  the  gases  surrounding  a  space  vehicle  as  it  enters  a  planetary 
atmosphere  have  lead  to  interest  in  the  subject  of  radiative  gas  dynamics 
which  combines  radiative  transfer  ana  fluid  mechanics.; 

In  gas  dynamic  problems,  radiation  appears  as  a  frequency  integrated 
flux  in  the  energy  equation.  Therefore,  the  spectral  distribution  of  the 
radiated  energy  is,  in  principle,  not  required.  In  practice,  however,  no 
exact  method  of  avoid  ng  the  computation  of  the  spectral  distribution  has 
been  devised  except  in  the  limiting  cases  of  optically  thin  and  optically 
thick  media.  A  widely  used  approximation  which  avoids  spectral  difficul¬ 
ties  is  the  "gray  gas"  approximation  which  treats  the  radiative  properties 
of  the  gas  (specifically  the  absorption  coefficient)  as  independent  of 
frequency.  The  gray  gas  has  been  the  model  for  much  of  tho  work  in  radia¬ 
tive  gas  dynamics  but  recently  some  progress  has  been  made  in  the  use  of 
nongray  absorption  coefficients. 

A  principal  difficulty  in  any  attempt  to  approximate  real  absorption 
coefficients  is  that  not  much  is  known  of  the  relative  roles  played  by  the 
various  radiating  mechanisms  even  in  simple  transfer  proLlems  and,  in 
particular,  the  Importance  of  spectral  lines  is  uncertain.  The  work  pre¬ 
sented  here  examines  the  radiative  processes  which  occur  in  gases  for  the 


X 


simplest  possible  case,  namely  a  gas  of  uniform  properties  in  local 
thermodynamic  equilibrium.  Since  it  happens  that  the  relative  importance 
of  the  various  radiative  processes  depends  strongly  on  the  etate  of  the 
gas  and  the  path  length,  this  investigation  was  carried  out  for  a  specific 
range  of  conditions  chosen  to  include  those  achieved  by  a  space  vehicle 
reentering  the  earth's  atmosphere. 

This  work  is  in  three  parts.  Part  1  describes  the  various  radiative 
processes,  discusses  how  the  corresponding  absorption  cross  sections  and 
occupation  numbers  may  be  calculated,  presents  approximate  methods  for 
computing  the  radiative  energy  transferred  by  spectral  lines  and  gives  an 
account  of  the  relevant  theories  of  line  broadening. 

Part  2  examines  some  characteristics  of  line  radiation  in  a 
bydrogenically  distributed  spectrum.  Some  fairly  simple  calculations  help 
to  reveal  the  factors  which  determine  the  distribution  of  relative  import¬ 
ance  among  the  stronger  lines.  Also  presented  are  new  methods  of  account¬ 
ing  for  the  many  weak  lines  and  an  indication  of  their  importance  relative 
to  the  strong  lines. 

Part  3  is  an  account  of  the  calculation  of  specific  intensities  in 
uniform  nitrogen  allowing  for  all  important  radiative  processes  and  with 
particular  emphasis  on  the  lines.  The  data  is  tabulated  and  discussed, 
the  method  of  calculation  is  described  and,  finally,  the  results  are 
presented  graphically  and  interpreted  in  detail. 
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i.  INTRODUCTION 

The  basic  equation  governing  the  transfer  of  radiant  energy  in  a 
non-scattering  gas  in  local  thermodynamic  equilibrium  may  be  written 

dl  /ds  s  K  (B  -  I  )  (1.1) 

V  V  V  V 

where  I  is  the  specific  intensity  of  radiation  in  the  direction  s  per 
unit  frequency  v  per  unit  time  per  unit  area,  is  the  Planck  distribu¬ 
tion  function  and  is  the  linear  absorption  coefficient.  The  quantity 

is  here  defined  to  include  stimulated  emission  and  completely  specifies 
the  radiative  properties  of  the  gas;  this  Part  is  concerned  with  the  methods 
of  calculating  Ky  which  are  currently  available  for  light  non-hydrogenic 
atoms  and  ions,  in  particular  nitrogen. 

1.1.  The  Absorption  Coefficient 

Absorption  and  emission  of  radiation  in  a  gas  correspond  to  changes  in 
internal  energy  undergone  by  atoms  and  molecules  and  changes  in  the  transla¬ 
tional  energy  of  free  electrons.  Since  the  discussion  of  this  Part  is  re- 
s.„  iCtea  to  atoms  and  ions,  the  energy  changes  we  have  to  consider  are 
exclusively  electronic.  Fortunately,  this  aspect  can  be  separated  from  the 
problem  of  finding  the  number  of  particles  in  a  given  state  so  that  the  ab¬ 
sorption  coefficient  due  to  a  species  S  can  be  written 

K  _  =  J  Nc.  o  .[1  -  exp( -hv/kT)]  (1.2) 

v ,s  .Sivi 
i 

where  N  .  Is  the  number  density  of  particles  of  species  S  in  state  i 

b  X 

(known  as  the  occupation  number  of  state  i),  o  .  is  tne  cross  section  for 

v  i 

absorption  from  state  i,  ana  the  factor  [1  -  exp(-hv/kT)]  accounts  for 
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stimulated  emission.  The  main  problem  in  finding  the  occupation  numbers  is 
to  determine  the  total  number  of  particles  present  for  each  species,  this  is 
considered  in  detail  in  Chap.  2  where  plotted  values  cf  the  species  composi¬ 
tion  of  nitrogen  are  presented. 

1.2.  Absorption  Cross  Sections 

The  cross  section  a  .  of  Eq.  (1.2) is  the  sum  of  a  number  of  cross 

vi  n 

sections  each  corresponding  to  an  independent  electronic  process  with  initial 
state  i  and  it  is  usually  most  convenient  to  consider  each  one  separately. 

A  sketch  of  a  simple  arrangement  of  electronic  energy  levels  is  shown  in 
Fig.  1.  As  can  be  seen  from  the  figure,  three  distinct  types  of  absorptive 
transition  can  occur:  bound-bound,  bound-free  (or  photoionization)  and 
free-free.  Free-free  transitions  give  rise  to  cross  sections  which  are  con¬ 
tinuous  for  all  values  of  frequency,  bound-free  cross  sections  are  zero  for 
frequencies  below  a  threshold  value  but  finite  and  continuous  for  higher 
values,  in  contrast  the  cross  section  of  a  bound-bound  transition  is  signifi¬ 
cant  only  over  a  very  small  frequency  interval  and  can  in  some  senses  be 
treated  as  a  singularity  occurring  at  a  single  frequency. 

Bound-bound  transitions  are  discussed  in  Chaps.  3  and  **:  Chap.  3  is 
concerned  with  the  calculation  of  the  frequency  integrated  line  absorption 
coefficient  while  Chap.  4  discusses  line  profiles  as  predicted  by  the  rele¬ 
vant  theories  of  line-broadening.  Finally,  Chap.  5  describes  methods  for 
calculating  the  cross  sections  due  to  bound-free  and  free-free  transitions. 
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2.  SPECIES  COMPOSITION  AND  OCCUPATION  NUMBERS 


2.1  Introduction 

In  a  gas  in  local  thermodynamic  equilibrium  the  occupation  number 
density  N^,  of  the  ith  energy  level  in  a  given  species,  S,  may  be 
found  from  the  Boltzmann  formula 

Ngi  =  g.Ng  exp(-Ei/kT)/(0)s  (2.1) 

where  is  the  total  number  of  particles  of  species  S  per  unit  volume, 

E^  is  the  energy  of  the  ith  level  above  the  species  ground  state,  g^  is 
the  statistical  weight  of  the  ith  state,  and  (Q)g  is  the  partition  func¬ 
tion  of  species  S.  The  species  composition  is  determined  by  a  set  of  mass 
action  equations  together  with  species  conservation  and  charge  neutrality 
conditions.  In  the  case  of  a  single  diatomic  gas,  like  nitrogen,  the  mass 
action  equations  are  ionization  and  dissociation  equations  which  have  the 
general  form 


(Q)S(Q)S,  _D/kT 


(2.2) 


For  an  ionization  reaction,  for  example.  S’  is  the  electron  gas  and  D  is 
X,  the  ionization  energy  of  species  SS'.  Since  the  partition  functions 
(Q)g,  (Q)g,  and  (Q)<.<.,  depend  on  the  number  densities  of  the  various 
species  (except  at  low  temperatures),  their  evaluation  is  coupled  to  the 
solution  of  the  mass  action  equations  (which  themselves  must  be  solved  by 
an  iterative  procedure).  It  is  clear  therefore  that  an  accurate  determina¬ 
tion  of  the  species  composition  is  a  complex  problem.  Section  2.2  discusses 
the  problem  of  predicting  the  truncation  point  of  a  partition  function. 
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Section  2.3  describes  the  calculation  of  species  exposition  in  the  particular 
case  of  nitrogen;  the  results  presented  here  will  be  of  use  in  Part  3. 


2.2.  Partition  Functions 

The  partition  function  of  species  S  may  be  written  to  a  close  approxi 

nation  as  the  product  of  a  translational  partition  function,  ^Qtr^s  s 
2  3/2 

V(2»mgkT/h  )  (where  V  is  the  volume  and  m^  the  mass  of  a  particle  of 
species  S),  with  an  internal  partition  function  (Q^  )  which  involves 

the  energy  levels  associated  with  the  internal  structure.  In  the  case  of  a 
diatomic  molecule  the  internal  energy  levels  are  those  of  rotational,  vibra¬ 
tional  and  electronic  excitation  whereas  only  the  electronic  levels  exist  for 
a  monatomic  particle.  Unfortunately,  the  higher  temperatures  at  which  the 
monatomic  particles  occur  lead  to  a  complication  of  the  electronic  partition 
function  which  is  discussed  in  the  next  few  paragraphs. 

The  definition  of  the  electronic  partition  function  of  species  S  is 


(Qel>S  = 


i=l 


gie 


-Ei/kT 


(2.3) 


where  g^  is  the  statistical  weight  of  the  ith  energy  level,  E^  is  its 
value  relative  to  the  ground  state  of  species  S  and  the  index  i  is  chosen 
such  that  E.  >  E..  It  is  well  known  that  the  series  in  Eq.  (2.3)  is 
divergent  since  it  represents  the  partition  function  of  an  atom  in  an  un¬ 
bounded  volume  of  non-interacting  particles.  When  the  interactions  between 

particles  are  allowed  for,  the  series  terminates  at  some  level,  i  , 

max 

Eq.  (2.3)  thus  becomes 
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(Q 


el^S 


i 

max 


l 


i=l 


-E./kT 

e 


(2.4) 


At  low  temperatures  (*  7000°K)  the  exponential  factors  in  Eq.  (2.4)  become 
so  small  that  the  summation  may  be  truncated  after  a  few  terms  and  the  value 
of  imax  not  important.  At  higher  temperatures,  on  the  other  hand,  later 
members  of  the  series  become  important  and  because  the  energy  levels  tena  to 
a  limit,  x  the  ionization  energy,  the  last  terms  in  the  series  are  frequent¬ 
ly  the  dominant  ones  since  the  g.  increase  with  i.  The  value  of  i  is 
then  of  great  importance  in  the  calculation  of  the  partition  function  and 
since  the  value  of  i  depends  on  the  interaction  force  and  the  density  of 
interacting  particles,  it  follows  that  the  partition  function  must  be  density 
dependent.  One  consequence  of  this  density  dependence  is  that  the  perfect 
gas  equation  of  state  does  not  strictly  hold  (since  the  derivative  of  Q 
with  respect  to  volume  is  not  zero)  but  it  is  shown  in  Appendix  I  that  the 
departure  over  the  range  of  conditions  considered  here  is  negligible  in  the 
case  of  hydrogen  (which  is  the  only  gas  which  can  be  treated  analytically). 

The  termination  of  the  partition  function  series  is  interpreted  as  a 
reduction  in  the  number  of  bound  states  to  a  finite  number  which  implies  a 
reduction  in  the  ionization  potential.  There  will  therefore  be  a  correspond¬ 
ing  change  in  the  value  of  x  to  be  used  in  the  mass  action  equation:  the 
vacuum  ionization  potential  x  has  to  be  replaced  by  an  effective  ionization 
potential  =  X  *  Ay .  This  effect  is  additional  to  the  influence  of  the 

partition  function  oi  the  species  composition.  The  effective  ionization 
energy  is  related  to  the  termination  of  the  partition  function  ,by  the  obvious 
relation  E.  <  X  while  for  hydrogen-like  states,  the  highest  bound 
quantum  number  is  given  by  the  largest  value  of  n  satisfying  the 


relation 
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hRc 


(n  ) 
max 


2  - 


<  Ax 


(2.5) 


where  R  is  the  Rydberg  wave  number  constant. 

The  depression  of  the  ionization  potential  clearly  depends  on  the  mech¬ 
anism  of  interaction  between  the  particles;  at  sufficiently  low  temperatures 
the  interaction  will  be  through  a  Van  der  Waals  force  but  at  such  temperatures 
(except  at  extremely  low  densities)  the  partition  function  will  be  dominated 
by  its  leading  terms  and  the  details  of  the  cut-off  point  are  not  important. 

At  higher  temperatures,  the  dominant  interaction  is  by  the  electrostatic 
forces  of  charged  particles  and  it  is  with  this  effect  that  we  are  concerned. 
The  problem  of  the  determination  of  Ax  as  produced  by  the  interaction  of 
charged  particles  has  received  considerable  attention  but  as  yet  there  does 
not  appear  to  be  complete  agreement  on  the  model  to  be  used.  According  to 
most  authors,  an  adaption  of  the  theory  due  to  Debye  and  HUckel  for  electro¬ 
lytes^-  is  valid  in  some  form  and  over  some  range  of  relatively  low  densities 

2 

and  high  temperatures  while  at  high  densities  Unsdld's  nearest  neighbor 
effect  is  dominant;  a  recent  paper  by  Stewart  and  Pyatt^  presents  a  more 
elaborate  theory  in  which  the  Debye-HUckel  and  nearest  neighlor  expressions 

,  4 

appear  as  limiting  cases.  Bond,  Watson  and  Welch,  however,  suggest  that  at 
low  densities  and  a  low  degree  of  ionization  the  Bohr  radius  should  no: 
exceed  the  mean  distance  between  particles.  These  theories  are  discussed  in 
the  next  paragraph  but  it  should  be  remarked  that  even  the  best  of  them  are 
not  highly  accurate  since  they  all  postulate  a  sharp  cut-off  between  bound 
and  free  states  and  neglect  the  displacement  of  the  higher  bcur.J  levels 
which  occurs  due  to  external  force  fields. 

The  Debye-HUckel  theory  examines  the  screening  effect  <-f  cha-^ed 
particles  surrounding  an  atom.  According  tc  this  theory  the  characteris* ic 


distance  for  the  reduction  in  the  Coulomb  potential  it  the  Debye  radius, 
Rp,  defined  by 


R  * 

kd  - 


kT 


4»e2(N#  ♦  l  t2Ng) 


1/2 


(2.6) 


where  N  is  the  number  of  electrons  per  unit  volume  and  N  the  number  of 
e  z 

z  times  charged  particles  per  unit  volume.  One  can  now  either^  calculate 
the  polarization  energies  of  an  ion  and  an  electron  and  fro/n  them  the  reduc¬ 
tion  in  ionization  energy  or^  treat  all  states  whose  semi-major  axis  exceeds 
the  Debye  radius  as  free.  These  two  procedures  are  qualitatively  equivalent 
but  lead  to  results  which  differ  by  a  factor  of  two.  These  are  respectively 

Axz  *  (z+De2/^  (2.7) 

and 

Axz  *  (z>l)e2/2RD  .  (2.8) 

The  Debye-HUckel  theory  involves  an  averaging  over  charges  which  is  only 
valid  if  the  Debyo  sphere  contains  a  minimum  number  of  charged  particles, 

7 

according  to  Duclos  and  Cam be 1  the  condition  is 


N 

e 


(2.9) 


which  implies  that  there  must  be  at  least  one  sixth  of  a  charged  particle  in 
the  Debye  sphere.  Combining  Eqs.  (2.6)  and  (2.9)  gives 


M  ♦[  i2H 
e  L  t 


*  ♦  7  N 

e  *•  z 
z 


{»,  ♦  l  *\)U1 

z 


i  ri  (")  J 

e* 


1/2 


(2.10) 
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and  since  the  ratio  (N  ♦  J  z2N  )/(N  t  [  N  )  <  2.6  for  kT  <  3 eV  and 

e  *“  z  e  L  z 
-6  z  z 

p/p  >  10  ,  then  Eq.  (2.10)  must  become 

o 


n  ♦  y  z2n 

e  L  z 
z 


< 


_1_ 

7* 


(2.11) 


It  can  thus  be  seen  that  the  equality  of  Eq.  (2.9)  represents  an  upper  bound 
on  charged  particle  density  since  the  Debye  volume  increases  more  rapidly 
than  the  charged  particle  density  and  hence  the  number  of  particles  in  the 
Debye  volume  must  decrease  as  the  number  in  unit  volume  in  space  increases. 

Q 

Ecker  and  Kroll  agree  with  the  Duclos  and  Cambel  result  except  for  a  factor 

4 

of  3/**  in  the  limiting  density.  Bond,  Watson  *nd  Welch  on  the  other  hand 

suggest  that  the  Debye  criterion  should  be  used  provided  R  >  r  ,  the  mean 

distance  between  particles  (not  electrons),  otherwise  the  ionization  limit 

should  be  depressed  to  the  first  state  whose  Bohr  radius,  a  ,  is  less  than 

n 

r  .  Since  we  have  that 

m 


(2.12) 


where  is  the  number  density  of  atoms,  the  condition  for  validity  of  the 

Debye  theory  becomes,  inserting  Eq.  (2.12)  into  the  relation  R  >  r  , 

D  -  w 


n  ♦  y  n 

a  *-  z 
z 


(2.13) 


In  contrast  to  Eq.  (2.9)  this  constitutes  in  upper  bound  on  the  number 

of  particles  in  the  Debye  sphere.  As  an  example  suppose  that  the  first  ior 

% 

dominates,  then  since  N  *  N+  and  all  ether  species  are  negligible,  the 
combination  of  Eqs.  (2.13)  and  (2.6)  with  (2.10)  gives 
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288* 
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N 

e 


<  i  (*£!)' 

-  *  \2J 


(2.m) 


which  is  a  very  small  range  of  conditions.  Most  recent  workers,  * 
however,  use  the  Debye  criterion  at  much  lower  densities  than  the  lower  limit 
of  Eq.  (2.14).  Me  have  to  choose  between  the  methods  and,  in  the  absence  of 
any  direct  evidence,  we  follow  the  majority  and  do  not  use  the  mean  particle 
distance  criterion.  The  question  remains  as  to  whether  the  upper  bound  of 
Eq.  (2.10)  is  ever  reached  in  our  range  cf  conditions.  Some  numerical  checks 
confirm  that  for  temperatures  above  1/2  eV,  Eq .  (2.10)  is  satisfied  up  to 
densities  of  at  least  10  times  atmospheric. 

A  calculation  of  species  composition  can  therefore  be  set  up  as  follows. 
Starting  with  arbitrary  truncation  points  for  the  partition  functions,  the 
system  of  mass  actior  plus  conservation  equations  are  repeatedly  solved  and 
the  results  from  each  solution  are  used  to  calculate  the  depression  of  ioni¬ 
zation  potential  according  to  Debye-HUckei  theory  and  hence  to  correct  the 
truncation  points  of  the  partition  functions.  The  solution  of  the  mass 
action  and  conservation  equations  for  given  partition  functions  is  itself 
an  iterative  procedure.  A  rapidly  convergent  scheme  is  presented  in  Appen¬ 
dix  IV.  The  composition  calculation  for  nitrogen  (some  details  of  which  are 
discussed  in  the  next  section)  has  been  programmed  for  the  IBM  360  computer 
and  proves  tc  be  a  rapidly  convergent  computation  (50  points  take  less  than 
3  minutes  j . 


2.3.  The  Composition  of  a  Nitrogen  Plasma 

The  application  to  nitrogen  of  some  of  the  ideas  of  the  previous  section 
will  now  be  discussed.  The  results  of  this  section  will  fce  used  in  the 


radiative  transfer  calculations  of  Part  3. 
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The  lowering  of  the  ionization  potentials  of  any  of  the  components  can 
be  significant  but  the  higher  ions  have  such  high  excited  energy  levels  that 
the  truncation  point  of  tha  partition  functions  is  net  relevant  under  our 
conditions.  We  therefore  account  for  the  effect  in  N  and  N+  (details  of 
the  energy  levels  >  sed  are  given  in  Appendix  II)  and  include  terms  in  the 
higher  ions  which  contribute  more  than  1%  (see  Table  2.1  for  details);  the 
calculated  partition  functions  of  N  and  N+  are  plotted  in  figs.  2  and  3. 
Debye-Kuckel  theory  in  the  form  Eq.  (2.7)  was  used  since  the  criterion  is 
somewhat  less  arbitrary  than  that  of  £«,.  (2.8). 

It  should  be  mentioned  that  use  of  an  approximate  partition  function 
was  investigated  and  the  details  are  given  in  Appendix  III.  It  seems  worth¬ 
while  to  give  a  summary  here,  however.  If  one  truncates  the  partition  func¬ 
tions  after  the  first  few  terms,  the  errors  in  number  densities  of  N  and 
Kf  can  be  considerable,  but  it  turns  out  that  the  electron  density  and  the 
occupation  numbers  of  the  various  energy  levels  are  much  more  nearly  correct 
The  electron  density  is  little  affected  because  at  low  ionization,  where  it 
is  most  sensitive  to  errors  in  (Q  .).,,  either  the  temperature  is  low  and 
therefore  the  cut-off  point  is  unimportant,  or  the  density  is  high  and  the 
"exact",  density-dependent  truncation  point  is  at  low  energy  levels,  i.e,, 
close  to  the  arbitrary  cut-off  chosen  for  the  approximate  partition  function 
The  occupation  numbers  are  insensitive  to  partition  function  because  where 
the  truncaced  approximation  is  very  badly  in  error,  ionization  is  high  and 
the  partition  function  cancels  out  between  the  Saha  equation  and  the 
Boltzmann  formula  (Eq.  (2.1)),  whereas  at  low  ionization,  the  approximated 
partition  function  is  more  nearly  correct.  A  detailed  investigation  of  this 
approximation  (Appendix  III)  shows  that  .  rrors  in  occupation  number  of  about 
SO'fc  <-an  occur  at  the  higher  temperatures  and  densities  of  the  range 
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leV  <  kT  <  3eV,  10  6atmos  <  p/p  <1  atmos.  The  investigation  also 

o 

revealed  that  the  occupation  numbers  at  high  temperatures  and  low  ionization 
are  sersitive  to  the  cut-off  point  of  the  partition  function  and  hence  even 
our  density-dependent  truncation  point  will  not  yield  very  accurate  results 
under  these  conditions. 

Finally,  we  consider  the  partition  function  of  molecular  species. 
Molecular  species  only  exist  at  relatively  low  temperatures  and  hence  there 
is  no  problem  of  truncation  in  the  case  of  their  partition  functions.  How¬ 
ever,  vibrational  and  rotational  states  also  exist  and  are  coupled  to  each 
other  and  to  the  electronic  energy  levels.  As  demonstrated  in  Appendix  III, 
we  will  achieve  the  greatest  accuracy  if  the  partition  functions  used  in  the 
equations  of  chemical  equilibrium  are  the  same  as  those  used  in  the  Boltzmann 
formula.  Because  of  this  we  are  restricted  in  our  treatment  of  the  molecu¬ 
lar  partition  functions  of  nitrogen  to  one  which  is  consistent  with  the 
absorption  cross  section  material  of  Allen10  which  we  use  for  the  radiative 
properties  of  the  band-systems  (see  Part  3).  Allen's  data  is,  as  far  as  we 
are  concerned,  a  cross  section  per  particle  in  the  lower  electronic  state 
(superscript  i)  of  the  transition.  We  therefore  have  to  find  (N1)  ,  the 
number  of  molecules  in  the  lower  electronic  state,  with  reasonable  accuracy. 
This  may  be  done  by  first  solving  the  equations  of  chemical  equilibrium  for 
the  total  number  of  molecules,  Ng  ,  using  the  truncated  partition  functions 
ir  the  form 

qs 

“Ws  5  .i  exp 

1=1 

and  then  using  the  Boltzmann  relation  in  the  form 

(n1)s  *  <si)s(Qi)s(9v)sNs  • 
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The  rotational  and  vibrational  partition  functions  used  are  those  given  by 
Allen, 


and 


Q*  =  Cl  -  exp  (-hcwVkT )] 
v  e 


(Q1)  =  kT/hcB1 

r  e 


-1 


(2.15) 


(2.16) 


where  u1  and  B1  are  spectroscopic  constants  taken  from  Allen's  report. 
These  expressions  treat  the  molecule  as  a  harmonic  oscillator  and  rigid 
rotor  respectively.  For  the  first  excited  level  has  an  energy  of  6.2eV 

and  the  summation  for  (Q.  )M  may  be  truncated  after  the  first  term,  for 

lOt  Nj 

Nj  on  the  other  hand  3  excited  levels  must  be  included.  Details  are  to  be 
found  in  Table  2.1. 

The  partition  functions  of  Table  2.1  were  incorporated  in  a  computer 
program  to  find  the  particle  densities  in  nitrogen.  The  mass  action  equa¬ 
tions  for  given  partition  functions  were  solved  by  the  iterative  technique 
discussed  in  Appendix  IV.  The  results  as  functions  of  density  and  tempera¬ 
ture  are  presented  in  Figs.  4  to  7.  After  these  calculations  were  started, 
a  paper  by  Drellishak,  Aeschliman  and  Cambel^1  became  available;  they  also 
use  density  dependent  partition  functions  according  to  Debye-Hiickel  theory 
but  choose  the  cut-off  according  to  the  alternative  criterion  (i.e.,  bound 
states  have  a  semi-major  axis  less  than  the  Debye  radius)  this  gives  rise  to 
a  somewhat  higher  value  of  partition  function.  Bearing  in  mind  this  differ¬ 
ence  and  different  molecular  partition  functions,  the  results  of  Ref.  11 
are  in  reasonable  agreement  with  ours 


TABLE  2.1 


Expressions  used  for  the  partition  functions  of  nitrogen 


(QintJN2 


kt  »  i0M8.34[l-exp<-0.274/KT)]-1  4  18.56 


+ 


exp( -3.18/kT)  , 

l-exp( -0. 300/JcfT  ' 


(Q.  =  4.03  *  103  *  - -  *2——-—— 

inc  N?  l-exp(-0.292/kT) 

^el^H  density  dependent,  see  Appendix  II  for  energy  levels. 
^el^N+  density  dependent,  see  Appendix  II  for  energy  levels. 

(QelV+  =  6  +  12  exp(-7 ,102/kT ) 

^el^N***  =  1+9  exp(-8.340/kT) 

(Qei)N++++=  2+6  exp(-10.003/kT ) 
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3.  BOUND-BOUND  TRANSITIONS 

3,1.  Introduction 

Bound-bound  transitions  occur  between  discrete  energy  levels  of  the 
atom.  They  would  therefore  occur  at  a  precise  value  of  frequency  if  it  were 
not  for  certain  perturbations  of  the  energy  levels  due  primarily,  in  our 
case,  to  neighboring  particles.  These  perturbations  are  discussed  in  the 
next  chapter  under  the  title  "Line-Broadening"  and  it  is  sufficient  to 
observe  here  that  for  all  cases  of  interest  to  us  the  resulting  line  width 
is  very  small  in  the  sense  that  the  spectral  structure  of  a  line  is  on  a 
quite  different  scale  of  frequency  variation  than  those  of  either  the  Planck 
function  or  the  continuous  absorption  coefficient  arising  from  other 
processes. 

This  difference  in  scale  leads  to  two  important  simplifications  of  the 
transfer  problem.  In  the  first  place,  the  Planck  function  can  be  taken  as 
constant  across  a  line  and  for  some  line  profiles  analytic  expressions  for 
line  intensity  can  be  obtained  (Section  3.3  discusses  some  relevant  cases). 
In  the  second  place,  although  the  transfer  problem  is  highly  non-linear  in 
absorption  coefficient,  an  effective  separation  of  the  frequency  integrated 
intensity  into  line  and  continuum  contributions  occurs  if  the  state  of  the 
gas  is  uniform. 

This  second  simplification  can  be  easily  demonstrated  as  follows.  The 
formal  solution  of  the  equation  of  radiative  transfer  (Eq.  1.1)  in  a  uniform 
medium  with  cool,  transparent  walls  is 

f  dvB  [l-exp(-K  s)] 

(V  v 

o 


I 


(3.1) 
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Now,  let  K  ,  and  K  „  be  the  absorption  coefficients  due  to  lines  and 
vL  vC 

continuum  respectively.  Then,  K  =  K  ,  ♦  K  _  and,  adding 

v  vL  vC 

B^[exp(-Kv(,s )-exp(-KyCs )]  to  the  integrand  of  Eq.  (3.1),  we  get 


I 


J  dvBv[l-exp(-KvLs)]exp(-K^,s) 
o 


•f 


dvBy[l-exp(-KvCs)] 

o 


Finally,  if  we  recognize  that  the  term  l-exp(-K  s)  is  negligible  compared 

V  L 

to  the  second  term  everywhere  except  in  the  neighborhood  of  the  unperturbed 
line  frequencies  v^,  we  can  write 

1  *  l  hi  «P<-Kvics)  *  lC  (3'2) 


where  we  define 


and 


XLi 


B  . 
vi 


dv[l-exp(-KvLs)] 


dvB  [l-exp(-K  s)] 
v  r  vC 

o 


This  result  is  apparently  well  known.  *  The  interpretation  is  that  near 
a  line  the  line  and  continuum  coefficients  each  reduce  the  intensity  associa¬ 
ted  with  the  other  but  the  overall  effect  on  the  continuous  radiation  is 
negligible  because  the  frequency  intervals  affected  are  very  small  while  the 

reduction  in  each  line  gives  rise  to  a  factor  exp(-K  ._s). 

r  viC 

A  further  simplification  of  line  transfer  problems  arises  from  the  fact 
that  the  value  of  the  frequency  integrated  line  cross  section  is  independent 
of  the  line  profile.  Because  of  this  the  cross  section  for  a  transition  from 
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a  lower  state  J  to  an  upper  state  J'  can  be  written  in  the  form 


e 


(3.3) 


where  L(v)  represents  the  line  profile  defined  so  that  £*L(v)dv  =  1  and 
fjj,,  the  f-number,  is  defined  in  terms  of  the  dipole  matrix  element  of 
quantum  mechanics  and  thus  depends  on  the  wave  functions  of  the  two  states 
involved.  It  may  also  be  interpreted  as  a  correction  to  the  classical  solu¬ 
tion  for  the  energy  radiated  by  a  harmonic  oscillator  and  for  this  reason  it 
is  frequently  called  an  oscillator  strength.  The  determination  of  f  and 
of  L(v)  are  largely  separate  problems;  the  former  is  discussed  in 
Section  3.2,  but  a  complete  chapter  (Chap.  4)  is  devoted  to  a  discussion  of 
line  profiles.  The  final  tireesections  of  this  chapter  consider  solutions  to 
the  transfer  equation  for  lines  with  some  common  types  of  profiles  and  a 
treatment  of  the  problem  of  the  merging  of  closely  spaced  lines  as,  for 
example,  in  multiplets. 


3.2.  Oscillator  Strengths 

The  oscillator  strength  of  a  line  is  equal  to  a  constant  times  the 

frequency  integrated  cross  section  as  can  be  seen  from  Eq.  (3.3).  In 

Ruasell-Saunders  coupling  (which  holds  for  the  stronger  lines  of  most  light 

or  medium  elements),  the  radial  and  angular  contributions  to  the  matrix 

.  14 

element  may  be  separated  and  the  oscillator  strength  written  m  the  form 


9»2m  v  , , ,  . 

f  ,  - ?  f  (  *  )  f(  :  )02 

JJ*  3h  g.  v  /  V  ’ 

J 


(3.4) 


where  g.  is  the  statistical  weight  of  the  lower  state  and  the  remaining 
factors  are  as  follows.  ( )  is  the  relative  strength  of  a  line  within 


a  multiplet,  it  depends  on  inner  quantum  numbers  J  and  J'  and  on  total 

angular  momentum  quantum  numbers  L  and  L’.  Values  of  /  (.  )  have  been 

compiled  by  Allen. (  '. )  is  the  relative  strength  of  a  multiplet  within 

a  transition  array.  An  early  paper  by  Goldberg^  gives  tables  of  a  related 

quantity  from  which  j  (  )  can  be  obtained  with  the  aid  of  a  normalizing 
17 

factor,  fortunately  most  of  Goldberg's  values  have  been  normalized  and 

1 5  18 

tabulated  in  the  books  by  Allen  and  Aller.  In  order  to  take  account  of 

the  possible  parentages  of  number  of  equivalent  electrons,  a  further  paper 

19 

by  Menzel  and  Goldberg  is  sometimes  required  (where  Ref.  19  does  not  apply 

20 

see  Kelly  and  Armstrong  ).  These  tabulations  are  convenient  to  use  but 
unfortunately  they  do  not  cover  all  the  cases  of  interest  and  for  the  excep¬ 
tions  it  is  necessary  to  turn  to  the  more  general  expressions  obtained  by 

21  5  22 

Rohrlicn  (which  are  also  presented  by  Griem  )  and  by  Shore  and  Menzel. 

These  expressions  take  account  of  the  parentage  of  equivalent  electrons  and 

they  give  both  "  (  '  )  and  '  (  ...  )  in  terms  of  the  Racah  coefficients 

5  23-25 

which  are  tabulated  in  a  number  of  places,  ’  the  most  extensive  of 

24 

which  appears  to  be  that  of  Nikiforov.  The  final  quantity  appearing  in 
Eq.  (3.**),  a,  is  the  radial  matrix  element  for  the  jumping  electron.  It  is 
the  determination  of  this  quantity  which  is  the  main  problem  in  f-number 
calculations;  r">me  of  the  methods  currently  available  are  discussed  in  the 
following  paragraphs. 

2  6 

The  most  popular  general  method  is  that  due  to  Bates  and  Damgaard 

5 

(which  is  also  given  in  detail  in  Griem 's  book  ).  They  observed  that  the 
most  significant  contribution  to  the  radial  matrix  element  o  frequently 
comes  from  a  region  in  which  the  potential  has  almost  reached  its  asymptotic 
Coulomb  orm.  They  therefore  use  hydrogenic  wave  functions  but  replace  each 
principal  quantum  number  n  by  an  effective  quantum  number  n*  determined 
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from  the  experimentally  measured  value  of  tne  energy  level.  Comparison  of 
this  method  with  experimental  results  and  the  predictions  of  more  accurate 
theoretical  methods  shows  good  agreement  when  the  jumping  electron  is  out¬ 
side  closed  shells  or  in  an  excited  level,  but  may  be  orders  of  magnitude 
in  error  when  the  jumping  electron  is  one  of  a  number  of  equivalent 
electrons.  Griem^  has  evaluated  f-numbers  using  the  Bates-Damgaard  method 
for  many  of  the  lines  listed  in  Refs.  27  and  28.  A  check  of  his  values  for 
some  lines  of  NI  against  the  results  due  to  Kelly  (see  below)  of  more 
accurate  self-consistent  field  methods  shows  good  agreement  in  the  visible 
and  infrared,  but  an  underestimate  of  approximately  an  order  of  magnitude 
for  the  ultraviolet  lines.  This  discrepancy  is  due  to  the  fact  that  ultra¬ 
violet  lines  originate  from  particles  in  a  low  state  of  excitation  where  the 
Bates  and  Damgaard  method  is  least  accurate.  Another  defect  is  that  this 
method  cannot  predict  values  for  transitions  between  states  within  the  same 
shell  (which  often  give  rise  to  very  strong  lines).  Finally,  a  modifica¬ 
tion  to  the  Bates  and  Damgaard  method  which  takes  some  account  of  the  non- 

hydrogenic  form  of  the  wave  functions  may  be  found  in  the  paper  by  Burgess 
29 

and  Seaton  whose  primary  purpose  is  to  develop  a  similar  method  for 
bound-free  transitions.  This  modification  removes  the  worst  inaccuracies 
of  the  Bates  and  Damgaard  method  (naturally  at  the  expense  of  greater 
complication)  but  depends  on  experimental  knowledge  of  energy  levels  which 
is  unfortunately  not  always  available. 

Nonhydrogenic  wave-functions  are  used  in  the  class  of  methods  known  as 

self-consistent  field  approximations.  Of  these,  the  most  widely  used  for 

accurate  calculations  is  based  on  the  wave-functions  of  Hartree  and  Took. 

It  is,  however,  very  lengthy  and  has  only  been  applied  to  a  few  lines  of  a 

30 

few  elements.  Kelly  has  used  it  to  calculate  f -numbers  for  several 
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important  lines  in  the  ultraviolet  spectra  of  nitrogen  and  oxygen;  these 
calculations  can  be  expected  to  be  accurate  to  about  20%.  Also  due  to 

62  j 

Kelly  is  an  extensive  list  of  values  of  a*  for  nitrogen  and  oxygen  which 
he  calculated  using  a  simplification  due  to  Slater  of  the  Hartree-Fock 
method.  Thanks  to  these  values  one  can  calculate  the  oscillator  strengths 
for  any  of  the  important  lines  of  nitrogen  and  oxygen  to  within  a  factor  of 
about  2  and  can  reasonably  expect  the  average  accuracy  of  all  the  strong 
lines  to  be  rather  better  than  that. 


3.3.  The  Curve  of  Growth  for  a  Single  Line  with  Wing  Profile  b/Cv-v^)13 
It  was  pointed  out  in  Section  3.1  that  the  equation  of  radiative 
transfer  for  a  line  is  simplified  by  the  narrowness  of  the  line.  For  a 
single  line  Eq.  (3.1)  therefore  becomes 


I  =  B 


[l-exp(-K  s)]dv 
v 


(3.5) 


Performing  the  quadrature  in  Eq.  (3.5)  clearly  requires  the  specification  of 
the  function  K^(v).  However,  when  K^s  <<  1  for  all  v  the  expression 


recuces  to 


B  s  K  d\ 
v  j  v 

°  o 


(3.6) 


which  is  a  particular  form  of  the  well-known  optically  thin  limit;  it  shcul 
be  observed  tnat  this  expression  docs  r.ct  depend  cn  the  profile  (compare 
Eq.  (3.3)).  There  are  a  number  of  profiles  for  which  the  integration  in 
Eq .  (3.5)  can  be  performed  in  closed  form .  However,  we  will  be  concerr.ec 
here  with  the  asymptotic  limit  of  large  path-length  for  a  particular  class 


of  profiles. 


The  normalized  profile  L(v)  of  Eq.  (3.3)  which  we  will  consider  has 


the  form 


h(v )  'v 


a  >  1 


for  w*vc  >y  u •  The  results  to  be  derived  are  not  new,  at  least  expressions 

for  the  important  cases  of  a  =  2.0  and  a'  =  2.5  can  be  found,  but  no 

derivation  appeal's  to  be  available  and  the  general  form  does  not  appear  to 

have  been  obtained  before.  For  sufficiently  large  path-length,  the  line  will 

become  heavily  self-absorbed  near  the  line  center  and  the  intensity  in  this 

core  region  will  be  B  „  independent  of  the  cere  profile,  the  only  effective 

o 

part  of  the  absorption  coefficient  will  then  have  the  form 


K 

v 


K  dv  *  b/|v-v  ja 
I  V  ’o' 

0 


Defining  the  transformation  z  =  K  s  and  taking  appropriate  crr.re  of  the 

v 

limits  we  get  from  Eq.  (3.5) 


I  =  3  2  -  (fcs!  K  dv) 

v  a 
o 


v 


1/a  f  -(1+1/a).,  -z., 

I  z  (1-e  )dz 

L 


(  3 . 7  ) 


The  integral  over  z  of  Eq.  (3.7)  may  he  performed  bv  parts  to  yielc 


1  =  2—  ( hs  i  K  dv)x/,aB  [ad-e  **)z  i"  4  I  ♦  a 

a  jo  V  vo  >c 


f  - 1 /a  - z 


z  e 


.  *1  /  'l  >  \ 

.  J  v  3  .  o  ) 


where  the  first  terra  in  the  square  bracket  is  zero  for  a  >  C  ar.c  the  secor.r 
term  is  a  form  of  the  complete  gawva  fvnet  ion  for  a  >  i.  uc«at  ic-r.  (?-.»•} 
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I  =  2B  (bsf  K  dv)1/a  rfl  -  -}  . 

"  L  v  1  a' 


(3.9) 


It  ^ow  becomes  convenient  to  define  quantities  w  and  t  as  follows 


*  b1/(a'1)[r(i  -  i))*''*'11 


(3.10) 


and 


t  =  s 


e.  dv/w 
v 


(3.11) 


Then  Eq.  (3.9)  can  be  written 


I/2wB  =  t 
vo 


1/a 


(3.12) 


and  the  corresponding  form  of  ti.e  optically  thin  result  (Eq.  (3.6))  is 


I/2wB  =  t/2  . 

v 

o 


The  intersection  of  these  two  curves  occurs  at 


(3.13) 


t  *  2 


a/(a-l) 


( 3 .  !*• ) 


For  t  >  2a/^a  we  call  the  line  self-absorbed  (or,  sometimes,  optically 
thi  .k),  for  x  <  2a^a  ^  we  describe  the  line  as  not  self-absorbed  or 
optically  thin. 

The  quantity  w  has  the  dimensions  of  frequency  and  can  be  interpreted 
as  representing  the  width  of  the  absorption  coefficient  profile  (as  we  will 
chow,  it  is  pre^isi»)y  the  semi-half -width  of  a  dispersion  profile).  The 
dijcensiohiess  parameter  t  is  an  optical  depth  based  on  an  absorption  cccff; 

•s* 

f- 

cient  j  K^dv/w.  The  relationship  between  this  absorption  coefficient  and 
that  of  a  line  with  a  dispersion  profile  (a  =  ?)  is  snatched  on  Fig.  8a;  it 


22 

can  be  seen  that  in  this  case,  the  value  of  K  dv/w  is  if  times  the  peak 

'0  V 

absorption  coefficient  of  the  line.  The  quantity  I/B  has  the  dimensions 

V0 

of  frequency  and  will  be  called  the  "effective"  or  "equivalent"  width  of  the 

line.  From  this  definition  it  follows  that  the  energy  transferred  by  a  line 

is  «_qual  to  the  effective  width  of  the  line  times  the  value  of  the  Planck 

intensity  at  the  line  center.  When  the  line  is  self-absorbed,  so  that  the 

central  part  of  the  line  is  blackened  out,  the  effective  width  has  a  fixed 

relationship  to  the  intensity  profile  (specifically,  half  the  effective  width 

is  the  distance  from  the  line  center  at  which  the  spectral  intensity  has  a 

value  of  1  -  exp{-[F(l  -  — )]”a)  times  the  black-body  value).  Figure  8b 

d 

shows  the  equivalent  width  and  the  spectral  distribution  of  intensity  for  a 
typical  self-absorbed  line.  When  the  line  is  not  self-absorbed,  the  equivalent 
width  is  still  defined  as  before  but  it  is  no  longer  meaningful  as  a  representative 
width  of  the  intensity  profile  (see  Fig.  8c).  Finally,  the  ratio  of  effective 
width  to  line  width  (i.e.,  absorption  coefficient  half-width)  is  I/2wBvq.  This 
quantity  has  the  value  21//^a”1^  at  the  intersection  point  of  the  asymptotic 
relationships  Eqs.  (3.12)  and  (3,13). 

The  two  cases  of  most  interest  are  those  of  a  dispersion  profile,  where 
a  =  2,  and  a  quasi-statically  broadened  line,  where  a  =  2.5.  The  corresponding 
forms  of  Eqs.  (3.12),  (3.13)  and  (3. .14)  are  to  be  found  in  Refs.  31  and  32. 
Reference  31  also  obtains  the  complete  curve  of  grow:h  of  a  pure  dispersion 
profile  and  it  is  of  interest  to  compare  the  asymptotic  forms  to  the  exact  curve. 

A  dispersion  profile  has  the  form 


L(v) 


,  ,2  2 
(v-v  )  +w, . 
o  dis 


(3.15) 


where  w..  is  the  semi-half-width  (i.e.,  it  is  half  the  width  of  the  line 
dis 
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where  the  absorption  coefficient  has  half  its  maximum  value).  This  profile 

becomes  w(j£s/ff^v~u0)2  *n  t^ie  w^nRs  and,  using  T(l/2)  =  in  Eq.  (3.10), 

we  obtain  w  =  w,.  . 

dis 

31 

'’’he  exact  curve  plotted  from  values  given  by  Penner  is  shown  on  Fig.  9 
and  compared  with  the  asymptotic  resu.  Cq.  (3.12)  and  (3.13).  It  can  be 
seen  that  the  line  is  well  approximated  by  its  asymptotic  forms,  the  maximum 
error  being  about  10%.  It  will  be  seen  later  that  in  practice  lines  some¬ 
times  depart  from  the  true  dispersion  profile  in  the  core.  In  such  cases  the 
asymptotes  discussed  here  are  a  slightly  better  approximation  to  the  correct 
curve  of  growth  than  the  true  dispersion  curve  of  growth  is. 

3.4.  The  Curve  of  Growth  for  a  Group  of  Lines  with  Wing  Profiles  b/(v-v^)a 

As  one  passes  along  a  ray,  the  intensity  and  hence  the  effective  width 
of  a  line  increases  and,  in  many  cases,  eventually  interferes  with  neighbor¬ 
ing  lines.  This  section  discusses  the  effect  cf  this  interference  as  it 
occurs  in  a  group  of  lines  in  a  uniform  gas.  The  method  of  presentation  is 
as  follows.  First,  we  consider  a  multiplet  of  lines  with  dispersion  profiles 
and  give  a  simplified  model  of  the  curve  of  growth,  we  then  extend  this 
result  to  a  more  general  group  of  lines  and  finally  apply  the  same  ideas  to 
a  group  of  quasi-statically  broadened  lines. 

A  multiplet  is  the  group  of  all  transitions  between  two  terms.  The 
individual  lines  therefore  arise  from  the  splitting  of  energy  levels  by  the 
inner  quantum  number,  J.  This  splitting  is  so  small  that  the  Planck  func¬ 
tion  can  be  taken  constant  across  the  multiplet;  the  multiplet,  however, 
cannot  in  general  be  treated  as  a  single  degenerate  line  since  for  an  impor¬ 
tant  range  of  conditions  the  individual  lines  are  isolated  from  each  other. 

An  exact  treatment  of  the  behavior  of  a  multiplet  must  take  into  account  the 
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spacing  and  strengths  of  the  lines  and  is  therefore  both  complicated  and 
different  for  each  multiplet. 

We  first  discuss  three  parts  of  the  curve  of  growth  of  a  typical 
multiplet.  Consider  a  multiplet  containing  N  lines  with  dispersion  pro¬ 
files,  and  assume  that  the  lines  each  have  the  same  half-width  w  but 
different  strengths.  We  define  an  equivalent  width  for  the  multiplet, 

W  =  I/B^  where  I  is  the  total  integrated  intensity  and  B  is  the  Planck 
o  o 

intensity  at  the  center  frequency  of  the  multiplet.  The  equivalent  width 

ratio,  W  ,  is  defined  by  W  =  W/2w  =  I/2wB  .  Now,  at  low  densities  and 

o 

small  path-lengths  all  the  lines  are  optically  thin  and  we  may  write 


W  =  T  /2 
r  ,1  m 


(3.16) 


where,  for  a  path-length  s,  t  is  given  by 
*  m 


N 

?  U  V- 

1=1  'o 


(3.17) 


Usually,  there  will  be  some  optically  thicker  condition  at  which  all  the 
lines  are  self-absorbed  but  no  significant  overlapping  occurs;  in  this  case 
one  gets 


W  *  ✓  N  /r, 
r  ,2  m* 


(3.18) 


where 


r?  -  !  (f  Vv)^( !  r  ^  ■  !  rv>r 

i=l  J o  i=l  1 o  i=l 


(3.19) 


and  is  a  property  of  the  multiplet.  One  can  show  without  much  difficulty 
A 

that  N  <_  N.  Finally,  for  very  large  optical  depths  the  multiplet  width 
will  be  much  less  than  the  effective  line  widths  and  the  whole  group  will 
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behave  like  a  single  line: 


W  .  =  ST  (3.20) 

r,3  m 

Clearly  from  Eqs.  (3.18)  and  (3.20),  the  error  in  treating  a  multiplet  as  j 
single  line  under  conditions  where  the  lines  are  isolated  but  self-absorbed 
is  an  underestimate  by  a  factor  J? .  If  the  lines  are  very  closely  spaced 
the  group  can  behave  like  a  single  line  as  •soon  as  it  becomes  self-absorbed 
but  multiplets  of  low-lying  states  invariably  exhibit  a  considerable  range  of 
the  intermediate  growth  rate  at  conditions  of  interest  to  us.  We  will  there¬ 
fore  consider  Eqs.  (3.16),  (3.18)  and  (3.20)  as  holding  for  three  major 
regions  of  growth  and  discuss  the  transitions  between  them.  A  typical  curve 
of  growth  is  sketched  on  Fig.  10. 

The  intersection  between  the  curves  defined  by  Eqs.  (3.16)  and  (3.18) 
occurs  at 

t  =  4N*  (3.21) 

m 


This  intersection  point  corresponds  to  the  point  of  change  from  optically 
thin  to  optically  thick  behavior  of  an  isolated  line  whose  absorption  coeffi¬ 
cient  is 


K  dv 
v 


N 

7  K  dv/N* 
v 


Curves  (3.18)  and  (3.20)  do  not  intersect  but  curve  (3.20)  will  start  to 
describe  the  growth  somewhere  in  the  neighborhood  of  the  optical  path  length 
for  which  the  effective  line  width  of  the  multiplet  treated  as  a  single  line 
(W3)  equals  the  frequency  spread  of  the  multiplet,  similarly,  it  will  depart 
from  the  behavior  of  Eq.  (3.18)  approximately  when  the  effective  line  width 
of  a  single  line  equals  the  average  inter-line  spacing  of  the  multiplet. 
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The  behavior  between  these  two  points  may  vary  considerably,  however,  as 
may  be  seen  from  the  following  two  examples. 

Firstly,  suppose  that  the  lines  are  of  uniform  strength  and  evenly 
spaced.  Then  the  overlapping  of  adjacent  lines  will  occur  simultaneously, 
the  entire  center-portion  of  the  multiplet  will  be  blackened  out,  and  growth 
will  occur  only  by  virtue  of  the  wings  of  the  two  outermost  lines.  There 
will  therefore  be  a  region  of  relatively  slow  growth  until  the  effective 
width  of  the  inner  lines  exceeds  the  spread  of  the  multiplet  and  the  behavior 
of  Eq.  (3.20)  commences,  this  type  of  transition  is  shown  on  Fig.  10. 

As  a  second  example,  suppose  that  the  N  lines  occur  in  M  well 
separated  clusters  of  closely  spaced  lines  of  approximately  equal  strength. 
After  the  growth  as  isolated  lines,  the  line  centers  of  the  M  clusters  will 
blacken  out  first  and  the  M  clusters  will  each,  as  in  the  previous  example, 
have  a  relatively  slow  rate  of  growth  until  a  point  is  reached  where  the 
effective  line  widths  exceed  the  inter-line  spacing  within  each  cluster  and 
the  multiplet  then  starts  to  behave  like  M  isolated  lines.  Finally,  after 
another  transition  region  of  slow  growth,  the  multiplet  will  behave  like  a 
single  line.  The  case  where  M  =  2  is  shown  on  Fig.  10. 

A  number  of  multiplets  of  NI  were  examined  and  found  to  be  more 
closely  approximated  by  the  case  of  evenly  spaced  lines  than  by  the  case  of 
isolated  clusters.  However,  the  concept  of  isolated  clusters  is  useful 
since  we  »/ ill  see  in  Part  3  that  it  is  sometimes  possible  to  group  together 
weak,  adjacent  multiplets.  In  most  cases  these  do  not  ever  merge  to  a 
super-line  in  our  range  of  conditions  and  the  model  we  treat  is  that  of  N 
lines  distributed  in  M  non-merging  clusters. 

The  intersection  of  the  optically  thin  region  with  the  self-absorbed 

but  isolated  region  does  not  depend  on  M,  so  for  t  <  4N*  we  use  Eq .  (3L16 ). 

m  - 
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For  t  >  4N*’  we  treat  the  lines  as  isolated  (Eq.  (3.18))  ^  ’  W  =  D  where 
m  i 

D  is  the  sum  of  the  widths  of  all  the  clusters.  Then  for  D/N*Vm* > >  D 
the  transition  region  of  slow  growth  commences  and  the  equivalent  width  is 
approximately  constant  and  given  by 


W23  *  D  (3.22) 

From  here  on  (W^  s  W^/M*/^  >  D)  we  treat  the  system  as  M  isolated 
clusters  by  means  of  the  equation 


(3.23) 


where 


I  ({ xvi«W1/2/(  I  {' 

i=l  J  V1  i=l  l 


K  .dv) 
vi  1 


1/2 


(3.24) 


Lines  with  moderately  highly  excited  upper  states  are  subject  to 

quasi-static  broadening,  as  discussed  in  the  chapter  on  line  broadening 

(Chap.  4).  They  are  also  relatively  weak  and  at  least  partially  degenerate. 

The  ideas  of  this  section  can  be  applied  with  certain  modifications.  No 

permanently  isolated  clusters  exist  so  that  the  equivalent  of  /M5^  is  unity. 

2/5 

The  growth  rate  is  slower  (a  single  self-absorbed  line  grows  as  t  com- 

1/2 

pared  to  the  dispersion  rate  of  t  )  and  the  equations  to  be  used  are  as 

follows.  We  start  by  recalling  from  Section  3.3  the  definition  of  equivalent 
line  width 


w 


b2/3[r(3/5)]5/3 


(3.25) 


where  b  is  a  constant  obtained  from  the  asymptotic  form  of  the  normalized 

5/2 

Holtsmark  profile,  L(v)  ^  b/(v-vQ)  .  The  corresponding  optical  depth  for 
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a  group  of  lines  is 


N  f" 

»  l  K  ^dv/w 

i=l  *  V1 


(3.26) 


and  the  effective  masher  of  lines  is 


5  C  I  (j  K  iv)J/5/(  f '|  K  dv) 
i*l  J  V1  i=l  J  V1 


2/5.j5/3 


(3.27) 


Then,  if 


T  <  25/3N* 

m  - 


(3.28) 


we  use 


I  *  B  t  w  . 
v  m 
o 


(3.29) 


C  . 

However,  if  t  >  2  n1* ,  the  expression  to  be  used  depends  upon  the 

in 


effective  width 


W  *  2(H*)3/5  t2/5  w  . 

10 


If  W  <  D  we  have 


I  *  B  W 

v 

o 


*  -3/5 

but  if  W  >  D  and  W(N  )  <  D  the  intensity  is  given  by 


I  *  B  D 
v 

o 


and,  finally,  for  W(N*)  3/^5  >  D  we  have 


(3.30) 


(3.31) 


1*3  W(nV3/5 

v 

o 


(3.32) 
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These  growth  expressions  have  been  programmed  as  a  FORTRAN  subroutine 
for  a  general  wing  decay  law  b/(v-vo)a  and  will  be  applied  in  the 
calculations  presented  in  Part  3. 


3.5  The  Curve  of  Growth  for  Doppler  Broadened  Lines 

This  section  is  concerned  with  the  radiation  from  Doppler  broadened 
lines.  The  approach  is  analogous  to  that  of  Sections  3.3  and  3.4.  However, 
the  form  of  the  Doppler  profile  causes  an  additional  difficulty  which  re¬ 
quires  special  treatment.  First  we  discuss  approximations  to  the  single- 
line  curve  of  growth  and  then  consider  the  treatment  of  groups  of  lines. 

The  broadening  of  a  line  due  to  the  thermal  motion  of  the  radiating 
particles  is  known  as  Doppler  broadening  and  is  discussed  in  Section  4.2. 

The  resulting  absorption  coefficient  can  be  written  in  the  form 

(v-v  )  /in  2  2 

Ky  =  Kv  exp{-  [ - - - 3  }  (3.33) 

o  u 


where  w^  is  the  semi-half-width  of  a  Doppler  broadened  line  given  by  (see 
Eq.  (4.2)) 


_  \  ,21 

c  'm. 


in  2) 


(3.34) 


and  Ky  is  the  peak  absorption  coefficient.  This  peak  value  is  related  to 
o 

the  frequency- integrated  absorption  coefficient  by  the  expression 


K  *  f  K  dv  *  /  hlLl 

V  I  V  w 

o  *o 


(3.35) 


The  intensity  is  given  by  substituting  the  expression  for  into  Eq. 

(3.5)  and  changing  the  variable  of  integration  to  x  =  (v-v^)  /in  2/w^,  to 
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obtain 


I  /Zn2 

B  w_ 
v  D 
o 


(1  -  «xp[  -  a  exp(-x  )]}dx. 
o  o 


(3.36) 


Tha  abova  intagration  cannot  ba  performed  in  closed  form.  Penner 
presents  a  series  representation  of  Eq.  (3.36);  the  series  is  uniformly  valid 
but  only  converges  reasonable  rapidly  at  comparatively  low  values  of  optical 
depth.  Penner  also  gives  a  simple  expression  which  is  approximately  correct 
for  snail  optical  depth  and  an  asymptotic  expansion  for  large  optical  depth. 
For  our  purposes,  the  approximate  expression  and  the  leading  term  of  the 
asymptotic  expansion  are  all  that  are  required.  These  are  respectively. 


^ 'd  •xpf-'V’J 

v  D 
o 


(3.37) 


R  w_  D 

voD 


(3.38) 


where 

*  K  s. 

D  -  v 

o 

1/2 

tq  is  analogous  to  the  t  of  Section  3.3  except  for  a  constant,  (In  2/w) 
(compare  Eq.  (3.35)). 

The  curves  defined  by  Eqs.  (3.37)  and  (3.38)  int«rsect  at  the  point 
Tp  *  26,  If  one  uses  Eq.  (3,38)  for  >  26  and  Eq.  (3,37)  for  <  26 
than  the  error  compared  to  the  exact  curve  is  always  within  about  15%.  The 
exact  and  approximate  curves  are  compared  on  Fig.  11. 


Using  Eqs.  (3.37)  and  (3. 38),  the  specific  intensity  from  a  multiplet 


of  N  lines  can  be  written  in  the  form 
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where 
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and  R^  and  are  given  by  the  expressions 
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By  reference  to  Section  3.4  it  can  he  seen  that,  for  a  dispersion  profile, 

/7 


respect- 


the  Quantities  corresponding  to  F.  and  are  unity  and 

X  z 

ivslv,  N  is  a  property  of  the  multiplet.  It  is  unfortunate  that  in  the 


Doppler  broadened  case,  both  R  and  depend  on  the  distance  s  and 

17 

on  the  line  width  w^,. 

Any  expressions  for  R^  and  would  be  useful  which  involved  only 

a  universal  property  of  the  multiplet  and  the  optical  depth  for  the  ■sultiniet 


T-.  .  Such  expressions  mean  that  it  isunnecassarv  to  treat  each  l*j*  lnriividu*liv . 

UJ* 
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Unfortunately,  it  has  not  proved  possible  to  express  R^  and  in  this 

way  except  where  the  lines  are  of  equal  strength.  In  this  case,  cne  gets 

R,  «  1  and  R_  *  N  ln(t.  /N)/tn  t_  .  (3.43) 

1  2  V  ,m  u  ,tn 

One  may  also  observe  that  in  the  general  expression  (3.41)  as  r  -*•  0, 

U  ,IH 

R^  ♦  1,  the  equal  strength  value  of  Eq.  (3.43), 

In  order  to  avoid  considering  the  individual  strength  of  each  line  of 
a  Doppler  broadened  multiplet,  Eq.  (3.43)  may  be  used  as  an  approx imat ion  for 
lines  of  unequal  strengths.  Some  checks  of  the  accuracy  of  Eq.  (3.43) 
against  actual  multiplets  revealed  errors  of  ?0t-30%. 

The  effect  of  overlapping  intensity  profiles  can  be  treated  in  exactly 
the  same  way  as  when  the  lines  are  dispersion  or  quasi-statically  broadened 
(see  Section  3.4). 

We  remark  in  closing  this  nection  that  the  far  wings  of  a  real  line 
never  have  a  Doppler  form.  However,  the  treatment  of  a  line  with  a  orofile 
formed  by  different  broadening  mechanisms  depends  on  the  relative  sizes  of 
the  associated  line  widths  and  a  discussion  is  delayed  until  Section  4.7. 
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4.  LINE  BROADENING 

4.1  Introduction 

It  has  been  pointed  out  in  Chap.  3  that  the  radiation  associated  with 
bound-bound  transitions  is  confined  to  a  small  interval  of  frequency  but 
that  when  self-absorption  occurs,  the  profile  L(v)  plays  an  important  role 
in  determining  the  amount  of  energy  transferred.  This  chapter  considers  the 
problem  of  calculating  the  function  L(v).  It  particularly  deals  with  forms 
which  are  simple  enough  to  be  of  use  in  the  prediction  of  radiative  transfer 
from  the  many  lines  of  a  non-hydrogenic  gas.  As  in  other  sections,  some  of 
the  discussion  is  about  nitrogen  and  will  be  applied  in  Part  3. 

There  are  three  mechanisms  which  can  cause  line  broadening,  (a)  Sir.'?'1 
atomic  states  have  finite  lifetimes,  the  uncertainty  principle  implies  that 
there  must  be  a  range  of  possible  energies  associated  with  each  state  and 
hence  a  range  of  frequencies  for  every  line  ■  phenomenon  is  called 
natural  broadening  (or  radiation  dampiry,*  (b)  Because  of  the  Doppler  effect, 
the  frequency  of  emission  (or  absorption)  of  radiation  by  a  particle  in 
motion  is  shifted  according  to  the  velocity  along  the  line  of  emission  (or 
absorption).  The  macroscopic  result  in  a  gas  where  the  particles  are  in 
thermal  motion  is  therefore  a  broadening  of  the  line,  (c)  The  third  mech¬ 
anism  is  a  perturbation  of  the  energy  levels  of  the  radiating  atom  resulting 
from,  interactions  with  other  part  icles ;  this  is  called  pressure  broadening. 

The  interactions  of  importance  are  due  to  the  van  der  Waals  force,  the 
electric  fields  of  charged  particles  (called  Stark  broadening)  and  the  coup¬ 
ling  of  two  similar  particles  (called  resonance  or  self-broadening) .  Of 
i  °  " 

We  are  indebted  to  Professor  Hans  R.  Criem  of  the  University  of  *»arvland 
who  very  kindly  discussed  with  us  a  number  of  points  in  connection  with 
the  treatment  of  line  profiles  and  who  supplied  us  with  references  37, 

3P  and  43. 
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these  meclianisms,  those  of  importance  in  a  non-hydrogenic  gas  in  the  range  of 
conditions  10  *  to  1  atmosphere  density  and  5000  to  35000°K  are  Doppler, 
resonance  and  Stark  broadening,  they  are  each  discussed  in  greater  detail 
in  the  following  sections,  Stark  broadening  is  given  the  greatest  attention 
because  of  its  dominance  and  complexity. 


4.2,  Doppler  Broadening 

The  theory  of  Doppler  broadening  is  well  established  and  may  be  found 
in  a  number  of  texts  ve.g..  Refs.  31,33).  In  the  absence  of  any  other  broadening 
mechanism  it  involves  a  statistical  average  over  the  Doppler  shifts  of  the 
individual  particles.  It  is  therefore  a  purely  temperature  dependent  effect 
which  produces  a  characteristic,  symmetrical,  narrow  line  with  very  low 
wings*  Analytically,  the  li^e  shape  is 


L(v)  =  — 

ft 


1  /in  2 


exp[~ 


(v-v  )2ln  2 


-3 


(4.1) 


where  vq  is  the  line  cento '  frequency  and  Wj, 
the  peak  intensity  (i.e.  L(vq+w^)  =  j  L(vq)) 


'll  /2VT 

WD  -  c 


in  2 


is  half  the  width  at  half 
it  is  given  by 


(4.2) 


where  m_  is  the  mass  of  the  radiating  particle. 


4.3  Resonance  Broadening 

Resonance  broadening  occurs  where  the  radiating  and  perturbing  particles 
ore  of  the  same  type  and  are  uncharged  (since  otherwise  Stark  broadening 
would  far  exceed  the  resonance  interaction).  From  the  point  of  view  of 
classical  mechanics  one  can  consider  the  gas  as  consisting  of  a  collection 
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of  oscillators  which  are  coupled  by  some  force  (e.g.  an  electrostatic  force), 
the  interaction  of  the  oscillators  through  this  force  results  in  each  oscil¬ 
lator  having  a  spiead  of  frequencies  around  its  natural  value.  When  the 
oscillators  have  the  same  natural  frequency,  a  resonance  effect  greatly 
strengthens  the  coupling  force  and  produces  a  greater  spreading  of  frequen¬ 
cies;  it  is  this  phenomenon  which  is  called  resonance  broadening.  Now, 
imagine  the  perturber  in  an  initial  state  1.  If  the  upper  state  2  of  the 
radiating  transition  can  be  reached  from  state  1  by  a  dipole  transition, 
then  the  perturbing  atom  will  broaden  state  2  by  a  resonance  interaction. 
Clearly  this  is  only  of  importance  where  there  are  a  significant  number  of 
particles  in  state  1.  This  implies  that  state  1  must  have  a  low-lying 
energy  level  because  there  will  only  be  a  significant  number  of  atoms  if 
the  temperature  is  low  while  at  low  temperatures  the  number  of  excited  states 
is  small. 

Resonance  broadening  results  in  unshifted  dispersion  profiles. 


L(v) 


1_ _ Wres 

V  (v-v  )2  +  w2 

o  res 


(4.3) 


where  w  (the  semi-half-width)  may  be  calculated  in  an  impact  approx ima- 
r*6s 
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tion  (see  next  section)  and  has  the  form 
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Here  f  ,  v  ,  g  and  g  all  belong  to  the  resonance  lire  and  are, 
respectively,  the  absorption  oscillator  strength,  the  frequency,  the  lower 
state  statistical  weight  and  the  upper  state  statistical  weight,  N  is  the 
number  of  perturbing  atoms  per  unit  volume.  Evaluating  the  constant, 
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£q.  (4.4)  becomes 
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where  is  the  Loschmidt  number. 

It  can  easily  be  checked  that  for  transitions  from  low-lying  states 
under  conditions  of  low  ionization  this  line  width  can  be  greater  than  that 
due  to  Stark  broadening  (see  next  section). 


4.4  A  Brief  Review  of  the  Theory  of  Stark  Broadening 

The  problem  of  Stark  broadening  (i.e.  the  broadening  of  spectral  lines 
by  charged  particles)  has  received  considerable  attention,  but  until  recently 
the  treatments  were  unsatisfactory  in  a  number  of  important  respects.  Nov, 
however,  the  work  of  Griem,  Baiunger,  Kolb  and  their  co-workers  has  pro¬ 
vided  a  general,  accurate  theory  (see  Refs.  5,  34  and  papers  referred  to 
therein).  This  section  gives  a  brief  description  of  some  of  the  main  fea¬ 
tures  of  the  broadening  mechanism  and  its  theoretical  treatment  with 
particular  emphasis  on  non-hydrogen ic  gases. 

A  detailed  discussion  of  the  Stark  effect  is  given,  for  example,  by 
35 

Bethe  and  Salpeter  and  so  only  a  very  brief  description  need  be  given  here. 
When  an  atom  is  subjected  to  an  electric  field,  the  energy  levels  are  each 
split  into  a  number  of  components.  If  this  is  viewed  as  a  perturbation  prob¬ 
lem  where  the  expansion  parameter  is  the  field  strength,  then  it  turns  out 
that  the  first-order  effect  is  identically  zero  and  the  magnitude  of  the 
splitting  is  proportional  to  the  square  of  the  field  strength  ("quadratic 
Stark  effect'  )  except  for  very  strong  fields  (where  higher  order  effects 
become  important)  and  for  energy  levels  which  are  degenerate  with  respect 
to  the  orbital  angular  quantum  number,  t,  where  the  first-order  (linear) 
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effect  occurs  and  is  dominant.  Since  hydrogen  and  its  isoelectronic  sequence 
(He  II,  Li  III,  etc.)  and  the  highly  excited  states  of  more  complex  atoms 
are  degenerate  with  respect  to  l,  the  linear  Stark  effect  governs  the 
splitting  of  the  associated  energy  levels;  under  these  circumstances  it 
turns  out  that  (in  contrast  to  levels  split  by  the  quadratic  Stark  effect) 
broadening  by  the  slowly  moving  ions  is  important  and  the  full  apparatus 
of  modern  line-broadening  theory  is  required  to  account  properly  for  the 
complex  line  shapes  which  result  (see  Refs.  5  and  34). 

We  turn  now  to  the  theoretical  treatment  of  line  broadening  through  the 
Stark  effect.  An  exact  analysis  would  proceed  as  follows:  the  field  at  the 
radiating  atom  is  calculated  for  a  general  sequence  of  perturbers  with  dif¬ 
ferent  *'elocities  and  trajectories,  this  field  is  then  substituted  into  a 
quantum  mechanical,  expression  describing  the  transition  and  a  statistical 
average  over  the  perturbers  performed  to  yield  the  line  shape.  Fortunately, 
in  practice  this  procedure  can  be  considerably  simplified  by  the  existence 
of  two  good  approximations  which  are  asymptotically  correct  in  the  opposing 
limits  of  interaction  times  which  are,  in  one  case,  very  much  less  than 
(|v-vq|)  (the  impact  approximation)  and,  in  the  other  case,  very  much 
greater  than  (lv“v0l)  *  (the  quasi-static  approximation).  In  the  impact 
approximation,  corresponding  to  rapidly  moving  particles,  the  process  is 
treated  as  a  series  of  discrete  encounters  with  single  particles.  In  the 
quasi-static  approximation,  for  slowly-moving  particles,  the  perturbations 
on  any  given  emitter  are  treated  as  coming  from  a  cloud  of  stationary  parti¬ 
cles.  In  both  cases,  statistical  averaging  is  necessary  to  achieve  a 
macroscopic  result.  Electrons,  being  rapidly  moving,  can  nearly  always  be 
treated  by  the  impact  approximation;  in  principle,  fcr  v-v^  sufficiently  large, 
electrons  should  be  treated  by  the  quasi-static  approximation.  Although 
this  effect  is  significant  for  some  hydrogenic  lines,  it  is  of  no  importance 
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for  non-hydrogenic  lines  whose  wings  decay  to  a  negligible  level  more  rapidly 
than  the  wings  of  hydrogenic  lines  (see  Ref.  5,  pages  92  and  93). 

Ions  sometimes  broaden  according  to  the  impact  approximation,  sometimes 
according  to  the  quasi-static  approximation  and  sometimes  according  to 
neither  limit.  Fortunately,  the  quasi-static  approximation  always  holds  for 
ion  perturbers  in  the  wingis  of  a  line  and  where  it  may  break  down,  i.e.  near 
the  line  center,  electron  broadening  usually  dominates.  Indeed,  the  influ¬ 
ence  of  ion  perturbers  in  a  non-hydrogenic  gas  may  often  be  neglected 
entirely  without  losing  reasonable  accuracy;  this  is  particularly  true  where 
the  radiating  particle  is  itself  an  ion,  since  the  local  density  of  positively 
charged  particles  will  be  reduced.  The  linear  Stark  effect  greatly  increases 
the  broadening  effect  of  ions  so  that  the  study  of  hydrogenic  lines  requires 
proper  accounting  for  both  electronic  and  ionic  broadening. 

The  simultaneous  action  of  ions  and  electrons  is  treated  by  a  simplified 
form  of  the  exact  treatment  described  earlier:  first,  the  Stark  splitting 
due  to  a  typical  field  of  stationary  ions  is  calculated,  then  the  electron 
broadening  of  these  lines  is  determined  from  the  impact  approximation  and, 
finally,  statistical  averaging  over  the  various  possible  field  strengths  is 
carried  out.  It  frequently  happens  that  in  the  upper  state  of  a  transition 
the  electron  is  so  much  more  weakly  bound  than  in  the  lower  state  that  pertur¬ 
bations  of  the  lower  energy  level  may  be  neglected.  In  general,  one  can  say 
that  this  is  the  case  for  all  transitions  except  those  where  the  jumping 
electron  is  one  of  several  equivalent  electrons  in  both  states  or  the  levels 
are  both  broadened  by  the  linear  Stark  effect. 

Finally,  a  few  words  relating  the  older  treatments  of  line  broadening  to 
modern  theory  will  prove  useful  in  the  next  section.  The  old  impact  (or 
collision)  theory  due  to  Lorentz,  Weisskopf,  Lindholm,  Foley  and  others  (see. 
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for  example,  Aller  for  a  description  and  references)  is  contained  as  a 
special  case  in  the  modern  impact  approximation.  The  older  treatment  is 
based  on  the  "adiabatic"  theory  which  treats  all  collisions  as  elastic  in 
the  sense  that  only  the  two  levels  between  which  the  system  is  radiating  are 
perturbed,  this  is  correct  only  for  low  perturber  velocities  (i.e.,  for 
electrons,  at  low  temperatures).  In  addition,  the  older  theory  does  not 
account  for  the  influence  of  overlapping  lines,  which  may  result  from  the 
near  degeneracy  of  the  field-free  spectrum  or  from  Stark  splitting  produced 
by  the  quasi-static  ion  field  or  even  from  the  excitation  by  the  ion  field 
of  nearby  forbidden  transitions.  The  current  form  of  the  quasi-static  approxi¬ 
mation  is  similar  to  the  old  statistical  (or  Holtsmark)  theory  except  that  it 
improves  the  distribution  function  for  the  positions  of  the  perturbing  par¬ 
ticles  by  accounting  for  their  mutual  interactions. 

4.5  Stark  Broadening  of  Strong,  Non-hydrogeric  Lines 

In  this  section  we  discuss  the  results  of  Stark  broadening  theory  as 
applied  to  strong  non-hydrogenic  lines  at  space  vehicle  re-entry  conditions 
and  with  particular  reference  to  nitrogen. 

According  to  both  the  adiabatic  and  the  general  impact  theories,  an 
isolated  line  has  a  dispersion  profile,  however,  the  shift  (d)  and  the  half¬ 
width  (w  )  are  different  in  the  two  theories.  We  recall  from  Section  3.3 
a  dispersion  profile  has  the  form 

w 

L(v)  =  7 - .  (4,6) 

(v-v  *d)2  *  w2 
o 

Before  discussing  expressions  for  w  ,  which  will  be  our  main  concern 
in  this  section  (shift  is  of  no  importance  in  radiative  transfer  problems), 
we  first  show  that  ion  broadening  may  be  neglected.  The  influence  of  ion 
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broadening  in  non-hydrogenic  gases  appears,  in  part,  in  an  asymmetric  line 
5  -2 

shape  such  that  one  wing  decays  as  (v-vq)  (as  in  the  absence  of  quasi- 
static  broadening)  while  the  other  wing  is  asymptotic  to  (v-vq)  •  Before 
this  latter  term  is  of  the  same  order  as  the  electron  contributed  impact 
broadening,  the  absorption  coefficient  of  a  typical  nitrogen  line  has 
dropped  to  lO”'*  times  its  peak  value,  as  may  be  checked  from  the  formulae 
and  tables  of  Griem  (ref.  5,  Chapter’  4).  The  influence  of  ion  broadening  on 
the  core  profile  can  also  be  estimated  from  Griem* s  work.  Equation  (4.90) 
of  ref.  5  is  an  empirical  expression  for  line  half-widths  due  to  both  elec¬ 
trons  and  ior.s  in  terms  of  the  half -width  due  to  electrons  alone  as  obtained 
from  the  inspection  of  calculated  helium  profiles.  This  equation  is 

wtotal  *  Cl  +  1*75a  (1  “  °‘75lflwei  (4.7) 

whrre  a  is  a  measure  of  ion  broadening  (absent  if  a  *  0)  and  r  is  a 
measure  of  the  mutual  interactions  between  the  perturbing  ions.  Taking  val¬ 
ues  of  a  from  Tables  4-5  and  4-6  of  Griem 's  book  and  estimating  values  of 
r,  one  finds  that  the  ion  broadening  effect  on  w  is  at  worst  (high  densi¬ 
ties)  about  10%  which  is  less  than  the  estimated  accuracy  of  the  best  theory. 

The  general  impact  approximation  results  in  an  analytic  expression  for 
line  width,  w,  (ref.  5,  Eqs.  4.68,  4.79  and  4.80)  which  is  accurate  to 

about  20%  but  which  unfortunately  requires  considerable  computation  to 

i - 

This  claim  requires  a  little  qualification.  It  does  not  hold  for  some  high 
orbital  angular  quantum  number  states  due,  as  discussed  in  the  next  section, 
to  Debye-shielding  of  the  electrons  and  partial  degeneracy  with  respect  to 
£.  although  the  experimental  results  of  Day  and  Griem39  for  two  4f  states  of 
Nil  show  that,  even  using  the  worst  estimate  of  experimental  errors,  the 
theoretical  results  of  ref.  5  are  correct  to  within  35%.  Another  source  of 
uncertainty  is  that  a  number  of  experimenters  (see  Jalufka,  Oertel  and 
Ofelt43  and  references  therein)  have  found  the  theory  to  be  in  error  by  a 
factor  of  2.5  to  3.0  for  All.  However,  the  results  of  Day  and  Griem39 
confirm  an  accuracy  of  20%  for  Nil. 
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evaluate.  In  view  of  the  other  inaccuracies  in  radiative  transfer  problems 
(in  oscillator  strengths,  energy  levels,  etc.)  such  computation  is  not  justi¬ 
fied  for  our  use  here.  However,  Griem^  has  evaluated  w  for  17  lines  from 
the  first  spectrum  of  nitrogen  and  75  lines  from  the  second  spectrum  of 
nitrogen  at  several  temperatures  covering  otir  range  of  interest;  these 
results,  normalized  by  their  values  at  20,000°K,  have  been  used  to  plot 
Tigs.  12  and  13.  The  temperature  dependence  of  the  line  widths  can  be  seen 
to  vary  considerably,  the  lines  with  negative  slopes  belong  to  high  orbital 
angular  momentum  upper  states  and  are  beginning  to  exhibit  the  dominance  of 
inelastic  collisions.  As  discussed  earlier,  Debye  shielding  of  the  electrons 
(not  accounted  for  in  Griem's  calculation)  may  be  significant  for  these 
states  and  therefore  the  line  widths  may  be  less  accurate  than  those  of  the 
lower  states.  Note  that  because  the  broadening  is  assumed  to  be  caused  en¬ 
tirely  by  perturbation  of  the  upper  states,  Griem's  results  provide  the 

widths  of  more  (unfortunately  not  many  more)  lines  than  those  actually 
>v 

tabulated.  Besides  choosing  an  average  temperature  dependence,  we  also 
have  to  obtain  an  expression  with  which  to  predict  the  line-widths  not  cov¬ 
ered  by  Griem's  calculations;  a  number  of  approximate  expressions  are 
available  in  the  literature  and  we  now  discuss  them. 

The  adiabatic  theory  of  line  broadening  due  to  upper  level  perturbation 

18 

gives  the  following  result  for  the  line  half-width 


e 


where  C4  is  the  quadratic  Stark  coefficient  of  the  upper  state.  In  its 
most  general  form,  is  given  by 

ince  the  influence  of  the  lower  state  is  neglected,  lines  with  the  same 
upper  state  should  have  the  same  frequency  width,  however,  the  value  for  Nil 
tabulated  in5  do  not  all  satisfy  this  requirement  —  apparently  due  to 
computational  errors. 


I 
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C4S 


e 


m  8ir2hc2 
e 


l  X2.f  . 
V  ui  ui 
1 


(4.9) 


or 

» 

Ctl  cm^sec-*  =  1.24  x  10  7  T  (X  .  cm)2  f  .  (4.10) 

4  jr  ui  ui 

where  the  summation  is  over  all  states  i  to  which  a  dipole  transition  •.both 
emissive  and  absorptive)  is  allowed  from  the  upper  state  u  whose  energy 
level  is  being  perturbed.  In  practice,  the  summation  is  usually  dominated 
by  a  few  strong  transitions.  Substituting  Eq.  (4.10)  into  Eq.  (4.8)  and 
evaluating  the  constant  gives 

weV:  1.3  X  I(f16  [|  l  (X  .  cm)2  f  .1  ]2/3  (N  cm_3)(kT  eV)1/6.  (4.11) 

?  UI  UI1  6 
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Griem  et  al  have  compared  Eq.  (4.11)  with  the  general  theory  for  several 

low-lying  lines  of  Hel  and,  over  our  temperature  range,  the  discrepancy  is 

within  a  factor  of  3.  It  will  also  be  recalled  that  the  adiabatic  theory  is 

asymptotically  correct  at  low  temperatures.  A  simplified  form  of  Eq.  (4.11) 

42 

is  available  for  a  hydrogenic  gas  since  then 

*  6 

2nC  =  l.u2  (—)  *  10~17  cm4sec"1  (4.12) 

Ct 


where  we  have  included  a  small  numerical  correction  pointed  out  by 
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Sibulkin  ;  n^  the  ffective  quantum  number  of  the  upper  state, 
stituting  this  re _  r\  Eq.  (4.8)  one  gets 


w  ( 


(N  cm'3)(kT  eV)1/6. 
6 


Sub- 


(4.13) 


No  derivation  of  Eq.  (4.12)  is  given  in  Ref.  42  but  presumably  it  is  obtained 
by  an  approximate  evaluation  of  Eq.  (4.10)  for  a  hydrogenic  spectrum  and 
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oscillator  strengths.  Since  this  treatment  cannot  account  for  same-shell 

transitions  it  is  likely  to  be  inaccurate  for  non-hydrogenic  gases. 

An  analytic  expression  is  also  available  in  the  opposite  limit  of  high 

34 

temperatures  where  weak,  inelastic  collisions  dominate.  Baranger  together 
with  Stewart  obtained  the  following  expression  for  the  width  of  non-hydrogenic 
ion  lines  due  to  weak,  inelastic  collisions  with  electrons. 


w  = 


3/3  m2 


N 


dv 


f(v>  I  _ J1 


.  T2T7TT  (4i>-1)oui«ff 


(4.14) 


where  the  summation  is  over  all  states  i  which  can  be  reached  from  the  per¬ 
turbed  state  u  by  dipole  transitions,  £  =  nvax(  l  ,£.)  and  g-,  is  the 
free-free  Gaur.t  factor  for  the  electron  transition.  If  we  now  put  g^  -  1, 

we  can  perform  the  velocity  integration  and,  assuming  a  Coulomb  force  law, 
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apply  a  sum  rule  to  obtain,  following  Armstrong  , 


2m  1/2  n 


i/3 


m 


N  (•)  (JL)  CSn  2  ♦  1  - 
2  e  'wkT'  v  z  '  u 


3£  (£  +1)] 
u  u 


(4.15) 


or,  numerically. 


eV  =  0.637  *  IQ'22  H  (kT  eV)‘1/2(— )  [5n*2  ♦  1  -  3C  (S  tl)].  (4.16) 

A  H  n  Ik 


U  u 


.46 


Using  an  approximation  due  to  Unsold  it  is  possible  to  reduce  Eq.  (4.16) 
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still  further  to  the  result  obtained  by  Stewart  and  Pyatt  , 


w  * 


( 2wm  ) 2 
e 


2wm  1/2  (n  )K 

W  \  — 


(4.17) 


cr 


*  4, 

»eV*  3.3  *  10'22  (kT  eV)‘i/2  (N  cm’3) 

e  ,2 


(4. IS) 
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& 

Note  that  Eq.  (4. IS)  contains  the  same  dependence  on  ny  as  does  Eg.  (4.13), 

the  z  and  temperature  dependences  are  quite  different ,  however. 
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From  the  work  of  Griem  et  al.  on  Hel  ,  we  expect  Eq.  (4.11)  to  agree 
with  the  more  accurate  results  of  Griem^  to  within  a  factor  of  2  or  3;  a 
number  of  direct  checteon  both  NI  and  Nil  confirm  this  expectation  and  we  can 
also  see  from  Figs.  12  and  13  that  a  temperature  dependence  of  7*^  is  an 
approximate  average  of  the  correct  behavior  (accurate  to  within  25%  fbr  our 
temperature  range).  Equation  (4.13)  agrees  with  Griem’s  results  to  withir. 
a  factor  of  5  for  NI  but  is  low  by  a  factor  of  approximately  40  for  Nil.  In 
contrast,  Eqs.  (4.16)  and  (4.18)  overestimate  the  widths  of  NI  lines  by  a 
factor  of  about  10  and  the  widths  of  Nil  lines  by  a  factor  of  4,  in  both 
cases,  the  temperature  dependence  of  Eqs.  (4.16)  and  (4.18)  is  incorrect  at 
theca  temperatures. 

For  low-lying  lines  not  covered  by  Griem’s  calculations,  the  best  that 
we  can  do  appears  to  be  to  use  the  result  of  adiabatic  theory,  Eq.  (4.11) 

(the  higher  lines  are  treated  quasi  -statically  as  discussed  in  the  next  sec¬ 
tion).  Thus,  we  use  the  following  expression  for  the  electron  impact  line 
width 

w  eV  *  A  N  /N.  (kT  eV)i/6  (4.19) 

e  L 

where  N^  is  the  Loschmidt  number  and  A  is  either  evaluated  by  patching 
Griem's  results  to  Eq.  (4.19)  at  20,C00°K  or  obtained  from  the  adiabatic 
expression  (see  Eq.  (4.11)) 

A  *  3.5  *  1C3  Ci  l  (A  .  cm)2  f  .j)2/3.  (4.2C) 

L,  ui  ui 

x 

There  is  one  exceptional  class  of  strong  lines  for  which  the  perturba¬ 


tion  of  the  lower  state  is  not  negligible  and  so  cannot  be  treated  bv  Eq. 
(4.20),  namely  same-shell  transitions  where  the  jumping  electron  is  one  of 
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several  equivalent  electrons  in  both  states  (which  normally  g i ve  lines  in 

1  g 

the  ultraviolet).  For  these  lines,  the  adiabatic  theory  Fives  (see  Aller^  ) 

A  =  3.5  *  103  [j  l  (X  .  cm)2  f  .  -  7  (X,.  cm)2  f  . |]2/3  (4.21) 

“  IU  U1  “  lx  U1 

1  s 

where  subscript  l  denotes  the  lower  state. 

4 . 6  Stark  Broadening  of  Lines  with  Highly  Excited  Upper  States 

Lines  with  highly  excited  upper  states  (not  covered  by  Hriem’s  results) 
will  be  subject  to  linear  Stark  broadening.  In  many  cases,  however.  they 
will  also  be  either  optically  thin  or  merged  into  a  pseudo-continuun,  so  tha 
the  energy  transferred  is  independent  of  the  profile.  However,  calculations 
carried  out  on  hydrogen  (see  Part  II)  suggest  that  under  some  important  con¬ 
ditions  (low  densities,  moderate  temperatures  and  high  path  lengths)  lines 
can  be  isolated  but  self-absorbed  up  to  upper  quantum  numbers  of  10  or  so. 

It  is  debatable  whether  the  accuracy  of  the  overall  calculation  would  be 
significantly  impaired  if  one  were  to  use  impact  line  profiles  for  these  nig 
lines  but,  taking  a  cautious  approach,  this  section  discusses  an  attempt  to 
account  for  the  peculiarities  of  these  lines. 

Perturbation  of  states  with  moderate  to  high  principal  quantum  numbers 
(n  >_  4  or  5)  ar.a  high  orbital  angular  quantum  numbers  is  complicated  bv 
certain  characteristics  of  these  states.  In  the  first  place,  L-S  coupling 

starts  to  break  down36  so  that  different  selection  rules  operate.  According 
rv  37  , 

to  way  ,  nowever,  this  should  have  little  effect  on  the  width.  Secondly, 
*r.e  energy  levels  become  relatively  cioseiv  spaced.  This  close  spac:",-;  has 
two  effects:  it  increases  the  importance  or  (..elastic  collisions  and  it 
implies  at  least  a  partial  degeneracy  in  orbital  angular  cuantu-t  number. 
Inelastic  collisions  are  due  to  electrons  whose  trajectories  are  relatively 
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distant  from  the  radiating  atom.  Such  electrons  are  subject  to  a  Debye 
38 

shielding  effect  .  This  shielding  is  not  taken  into  account  in  non- 

hydrogenic  electron  impact  theory  with  the  result  that  the  line  broadening 

is  overestimated.  However,  the  electron  impact  line  width  cannot  be  greater 

than  the  electron  impact  line  width  of  a  corresponding  line  in  a  hydrogenic 
38  39 

gas  *  and  we  thus  have  an  upper  bound  which  is  approached  as  the  degen¬ 
eracy  of  the  levels  increases.  A  further  consequence  of  the  degeneracy  of 

the  states  under  discussion  is  that  ion  broadening  becomes  important  in  the 

40 

line  wings.  Very  accurate  calculations  of  profiles  have  been  performed 
for  the  strongest  hydrogen  lines  and  approximations  valid  for  higher  hydro- 
gen  lines  are  also  available  .  Even  the  approximations  (obtained  by  folding 
a  dispersion  profile  due  to  electron  impacts  into  a  quasi-static  profile  due 
to  ion  fields)  are  too  complex  for  us  to  use  here.  Instead  we  take  account 
only  of  the  behavior  in  the  line  wings  (this  is  (see  next  section)  what  we 
do  for  the  dispersion  profiles  but  in  this  case  the  approximation  is  much 
worse).  Even  the  wing  shape  is  not  simple^  since  there  is  a  gradual  break¬ 
down  in  the  validity  of  the  impact  approximation  for  electrons  as  v-vq 
increases  until  eventually  they  too  are  described  by  the  quasi-static  theory. 

When  the  impact  approximation  holds,  the  line  wings  are  made  up  of  an  elrc- 

-2 

tron  contribution  which  decays  as  (V_VQ)  plus  an  ionic  contribution  which 

-5/2 

decays  as  (v-v  )  whereas  when  the  electrons  broaden  auasi-statically , 
o 

-5/2 

both  contributions  decay  as  (v-v  )  .  The  latter  is  the  simplest  to  deal 

o 

with  and  is  valid  at  the  lower  temperatures.  At  the  upper  temperatures  of 
our  range  the  intensity  profiles  of  the  lines  will  almost  always  merge 
before  the  effective  widths  reach  the  pure  quasi-static  region  but  nonetheless 
the  most  generally  applicable,  simple  treatment  appears  to  be  to  treat  both 
electron  and  ion  broadening  as  quasi-static.  The  analytic  expressions  for 
the  widths  of  degenerate  lines  are  obtained  as  described  in  the  next  paragraph. 
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In  principle,  in  the  theory  of  quasi-static  broadening  the  interaction 
between  the  perturbing  particles  should  be  taken  into  account.  This  effect 
is  small,  however,  for  neutral  radiators  except  at  high  densities  and  is  in 
any  case  least  in  the  wings  of  the  line  which  are  caused  by  close,  two- 
particle  encounters.  For  our  purposes  the  Holtsmark  theory  of  statistically 
independent  perturbers  will  be  sufficiently  accurate  even  for  charged  radia¬ 
tors.  The  simplest  form  of  this  theory  is  given  for  a  single  class  of 

18  42 

perturbers  by  Aller  and  by  Margenau  and  Lewis  .  The  resulting  asymptotic 

32 

form  may  be  written  to  conform  with  Griem  as 


L(v)  ■V  C(2irc  F  )3/2  —  - i-_  for  larf>e  Av(=v-v  )  (4.22) 

°  X*  <2„Av>5/2 

where  C  is  a  constant  for  a  given  line,  A  is  the  line  wavelength  and  F^ 
is  the  Holtsmark  normal  field  strength  given  by 

F  =  2.60  zeN2/3  (4.23) 

o 

where  N  is  the  —  .*v»er  of  perturbers  (whore  charge  must  bt  z)  per  cubic 

centimeter.  If  one  now  allows  for  more  than  one  class  of  perturber,  Eq. 

(4.23)  holds  but  F  becomes 

o 


F 

o 


2.60e  (H  +  I  z  N  )2/3 
6  z  P  ^ 

P 


(4.24) 


where  N  is  the  number-density  of  z  -times  ionized  particles.  Thus 
2P  P 

substituting  Eq.  (4.24)  into  Eq.  (4.22)  we  get 


L(v)  %  C(5.20ewc)3/2  (N  +  J  z  N  )  i - 

6  v  0  Zp  *3  (2wAv)5/2 


(4.25) 
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which  is  consistent  with  Griem 's  asymptotic  result  for  quasi -static  broadening 
by  electrons  and  singly  ionized  particles^.  The  constant  C  can  be  found 
from  the  approximate  analysis  due  to  Griem  for  Holtsmark  broadening  of 
hydrogen-like  lines.  If  one  follows  Griem' s  analysis  but  does  not  replace 
the  wavelength  by  a  hydrogen ic  value  one  gets 


_L  t  h 

16  1  „  2 
z  m  e4tr*c 
r  e 


■)3/J  >5  (n2  -  nf)3/S 

u  * 


(4.26) 


where  zp  is  the  charge  on  the  radiating  particle,  n^  is  the  upper  state 
quantum  number  and  is  the  lower  state  quantum  number.  This  result 

assumes  that  both  levels  are  degenerate  in  orbital  angular  quantum  number 
l  and  takes  account  of  the  perturbation  of  both  levels.  It  is  clear  from 
Griem' 8  analysis  that  where  the  lower  state  is  non-hydrogenic  (and  therefore 
may  be  ignored)  we  may  simply  set  n^  *  0  in  Eq.  (4.26).  Combining  Eqs. 
(4.26)  and  (4.25)  we  get 


L(v)  = 


_9_  f  1.30  h^ 

16  '  m  if  ' 
e 


-3/2 


(n2-n2)3/2 

(*e  + 1  \\  )  -~~  vr 

e  z  p  Zp  (2irAv  j  ' 2 
P 


(4.27) 


or,  numerically. 


L(v)  eV'1  =  2.12  x  10‘4  z  _3/2 

r 


,  <n2-n2)3/J 

(n  +  l  z  N  )  rp - - — - — =7^ 

6  z  P  Zp  L  (Ahv  eV)^ 

P 


(4.28) 


where  N.  is  the  Loschmidt  number.  The  corresponding  expression  for  the 
line  intensity  is  given  in  Section  3.3  where  a  width  parameter  is  introduced, 
namely 


5  b2/3  [r(3/5)]S/3 


(4.29) 
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5/2 

where  b  is  related  to  L(v)  by  L(v)  ^  b/(v-vQ)  .  Identifying  b  in 
Eq.  (4.28)  and  inserting  it  into  Eq.  (4.29)  yields 


w  =  *(*  /H  ♦  i  z  N  /N  )J/3  (4.30) 

ZP  P 

where 

A  5  0.691  x  lo"2  z  _1  (n2-n2).  (4,31) 

r  u  t 

« 

4.7  The  Relative  Importance  of  Doppler  and  Stark  Broadening 

In  this  section,  we  first  examine  the  curve  of  growth  of  intensity  of  a 
line  broadened  simultaneously  by  Doppler  and  any  dispersion  profile  produc¬ 
ing  effect.  The  important  parameter  is  the  ratio  of  dispersion  to  Doppler 
half-widths  which  is  next  calculated  for  some  typical  cases  in  nitrogen. 
Finally  from  these  results  it  is  demonstrated  that  Doppler  broadening  may  be 
neglected  over  a  wide  range  of  conditions. 

A  wj.aeiv  used,  approximate  expression  for  a  line  broadened  simultane¬ 
ously  by  the  Doppler  effect  and  a  dispersion  profile  producing  mechanism  is 


K  .  K  if’  £S±2ll  at  (4.32) 

v  0  ’  J-  a2  .  <5-t>2 


where 
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w  is  the  dispersion  semi-half-width  and  wD  is  the  Doppler  semi -half -width 
It  can  easily  be  shown  that 
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(4.33) 


which  is  the  pure  dispersion  profile.  The  first  limit,  a  ■+  “,  is  the  case 
of  dispersion  broadening  much  greater  than  Doppler  broadening  which,  of 
course,  must  be  asymptotic  to  the  pure  dispersion  case.  The  second  limit, 

C/a  ■*  •»,  shows  that  far  out  in  the  wings  of  the  line  C(v-vq)  »  w] 
the  absorption  coefficient  behaves  as  though  there  were  no  Doppler  broadening 
present;  this  result  is  to  be  expected  since  the  Doppler  profile  falls  off 
much  more  rapidly  than  the  dispersion  profile. 

Since  the  line  center  is  blackened  out  ut  large  optical  depths,  the 
detailed  shape  of  the  core  is  then  irrelevant  and  the  asymptotic  growth  is 
that  of  a  pure  dispersion  profile  with  semi-half-width  w,  and  (see  Sec¬ 
tion  3.3)  the  combined  intensity  is  given  by 


I  *  2B  /w  /"  K  d  s  for  large  s.  (4.34) 

v  0  v  v 

o 

On  the  other  hand  for  small  optical  depth,  the  intensity  is  independent  of 
the  line  shape  and  thus,  if  the  influence  of  Doppler  broadening  is  confined 
to  a  sufficiently  narrow  core  region,  the  curve  of  growth  will  be  independ¬ 
ent  of  the  Doppler  broadening.  The  curve  of  growth  of  a  combined  Doppler- 
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dispersion  profile  has  been  plotted  by  Penner  as  a  function  of  a  and 
optical  depth  based  on  the  Doppler  part  of  the  profile.  Figure  14  shows 
this  curve  of  growth  transformed  to  an  optical  depth  based  on  the  dispersion 
part  of  the  profile.  It  can  be  seen  that  for  a  >  0.5,  the  pure  dispersion 
growth  is  followed  almost  exactly.  For  a  >  0.05,  the  pure  dispersion 
curve  underestimates  the  combined  curve  by  at  most  a  factor  of  2  and  that 
only  for  a  small  range  of  optical  depths.  Values  of  a  which  we  may  expect 
in  our  problem  are  discussed  next. 
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iJsir.p  the  hydrogenic  approximation  to  the  adiabatic  result,  Eq.  (*+.13), 
ar.fi  the  Dnpoler  half-width  Eq.  (4,2)  we  get  for  a, 

a  =  26  *  (n")4  M  /N  (hv  eV)"1(kT  eV)“1/3.  (4.35) 

u  e  L  o 

Except  for  the  ultraviolet  lines,  the  important  lines  of  NT  lie  in  the 

it  . 

frequency  range  1  to  6eV  with  values  of  (n^)4  >  25,  hence,  for  kT  <  3eV 

a  >  75  N  /N.  . 
e  L 

.4  \  _4 

At  a  density  of  10  atmospheres.  N  /w  s  2  *  10  and  the  correspondir.r 

■  e  5 

-4 

value  of  a  is  0.015.  Thus  we  can  say  that  for  densities  greater  than  10 
atmospheres  the  maximum  error  on  any  line  is  a  factor  of  4  (with  the  excep¬ 
tion  of  some  ultraviolet  lines  where  the  factor  can  reach  8),  and  this  will 
only  apply  to  those  lines  (if  any)  whose  optical  depths  lie  in  the  range 
where  the  Doppler  profile  influences  the  radiated  intensity.  Except  in  the 
case  of  very  lonp  path-lengths,  most  lines  will  be  optically  thin  at  low 
densities  and  so  independent  of  the  profile  shape.  Even  the  errors  in  the 
ultraviolet  lines  are  not  as  severe  as  the  factor  of  8  obtained  above  since 
these  lines  are  subject  to  resonance  broadening  which  increases  the  disper¬ 
sion  half-width  above  the  value  used  in  Eq.  (4,35).  The  overall  errors  ir. 
neplectinp  Doppler  broadening  in  NI  are  therefore  likely  to  be  small 
compared  to  those  arising  from  other  sources. 

The  ionic  lines  will  suffer  more  on  account  of  their  higher  frequencies, 
indeed  the  resonance  lines  of  NUT  are  likely  to  be  Dopnler  broadened 
until  well  into  the  wings.  It  was  at  first  thought  that  the  lines  of  NTU 
would  not  carry  much  enerry  on  account  of  their  high  frequencies.  However, 
early  results  from  detailed  calculations  (reported  in  nart  3)  demonstrated 
that  a  number  of  these  lines  are  the  dominant  transfer  mechanisms  at  hirh 
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temperatures.  It  therefore  became  necessary  to  take  account  of  the  Donnler 
broadening  of  some  o*  the  ionic  lines.  We  do  this  as  rollovs. 

For  a  >  0.05  we  ipnore  the  Doppler  effect  and  treat  the  line  as  a  pure 
dispersion  line.  For  a  <  0,05  we  compute  the  intensity  due  to  a  pure 
Doppler  shape  by  the  approximate  methods  discussed  in  Section  3.5.  Next* 
we  calculate  the  intensity  accounting  only  for  the  dispersion  contour. 
Finally,  we  compare  the  two  values  of  intensity  and  take  the  larrer. 

The  treatment  just  described  is  equivalent  to  approximating  the  line 
by  a  pure  Doppler  core  region  followed  by  pure  dispersion  winps  and  nepiect- 
inp  the  transition  region  between  these  two  profiles. 

The  effect  on  intensity  of  such  an  approach  ir  shown  by  a  dashed  line 
on  Fig.  14  for  the  case  of  a  =  5  *  10*^.  It  can  be  seen  that  Lhe  error  in 
the  transition  repion  is  small. 
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5.  CONTINUUM  CROSS  SECTIONS 


5.1  Introduction 

This  chapter  discusses  some  of  the  methods  available  for  calculating  the 
cross  sections  for  bound-free  and  free-free  transitions;  in  the  calculations 
reported  in  Part  3  molecular  bands  are  included  in  the  continuum  but  since 
the  method  of  treatment  is  peculiar  to  the  band  systems  of  nitrogen,  the 
discussion  is  deferred  until  Part  3.  Continuous  cross  sections  are  more 
readily  available,  in  general,  than  those  for  lines,  and  the  associated 
radiation  has  been  more  extensively  studied.  This  means  that  the  task  of 
preparing  data  for  radiative  transfer  calculations  is  much  easier  for  the 
continuum  contributions  than  it  is  for  the  lines. 


5.2  Bound- Free  Cross  Sections 

The  bound-free  absorption  cross  section  for  a  particle  in  initial  state 
29 

i  may  be  written 


(  \  8  n£ 

°v  i  '  3  he 


«Vhv-S- 

°  gj 


(5.1) 


where  aQ  s  h2/4x2mc2  is  the  first  Bohr  radius,  is  the  statistical 

weight  of  the  initial  state  and  S  depends  on  the  transition  probabilities 
of  all  transitions  from  i  to  ionized  states  such  that  hv  1^,  where  I. 
is  the  ionization  energy  of  state  i.  In  Eq.  (5.1),  hv  is  measured  in 
Rydberg  units. 

In  Russell -Saunders  coupling,  the  radial  and  angular  contributions  may 
be  separated  and  the  expression 


(5.2) 
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is  obtained  where  the  primes  denote  the  final  states  of  the  electron  and  ion 
plus  electron  and  , ,  ,  is  the  integral  over  the  radial  component  of  the 
wave  function  (analogous  to  the  a2  of  Eq.  (3.4)).  Fortunately,;  ^ ^ ,  is 
generally  independent  of  L*  (Ref.  5,  page  109)  because  the  integral  is 
dominated  by  the  wave  functions  in  the  outer  regions  where  the  term  split¬ 
ting  is  small.  Multiplet  splitting  of  the  ground  state  may  also  be  ignored 
since  even  the  hydrogenic  decay  from  the  photoelectric  edge  (which  is  as 
1/v3  and  which  is  more  rapid  than  the  non-hydrogenic  decay)  is  responsible 
for  a  change  in  absorption  coefficient  of  only  about  1%  across  a  typical 
multiplet  width.  We  can  therefore  simplify  Eq.  (5.2)  to 

( 2Stl ) ( 2Ltl )  =  Cl-l^t-l  *  C£+lH.ltl  (5,3) 

where  C  ,  5  £  C  t .  When  the  jumping  electron  is  equivalent  in  the  ini- 
1  L' 

tial  state  to  q-1  others,  a  number  of  photo ionization  edges  will  normally 
result,  each  corresponding  to  a  different  ion  term  and  thus  to  a  different 
parentage  of  the  initial  state.  It  is  then  necessary  to  treat  each  parent¬ 
age  separately  and  apply  the  fractional  parentage  coefficients  mentioned  in 
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Section  3.2,  i.e.  one  gets 


c«-i  *  5HT  1  c.-l  *  <rsL(<'s'L')!2 


1  +  1 
2ftl 


(5.4) 


where  (Fc  (qS'L1))2  is  the  fractional  par. tag*  coefficient  to  be  found 

ub 


in  Refs.  19  and  20.  It  remains  therefore  to  calculate 


8  t ,  «ht((»\n,v,l)  (5.5) 

and  the  remainder  of  the  section  is  devoted  to  this  problem. 

As  in  the  case  of  line  cross  sections,  t>~e  most  accurate  method  of 
calculation  is  by  means  of  self-consistent  field  functions  but,  as  mentioned 
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before,  such  calculations  are  very  lengthy,  fortunately  thev  are  not  so 

necessary  in  the  case  of  bound-free  transitions  owing  to  the  widely  used 

and  fairly  accurate  method  of  Burgees  and  Seaton.  This  method  is  similar 

26 

in  approach  to  the  Bates  and  Damgaard  analysis  of  bound-bound  transitions. 
That  is,  it  is  a  semi-empirical  method  based  on  hydrogenic  wave  functions. 

It  is,  however,  more  complicated  to  apply  than  the  Bates -Damgaard  method 
because  it  includes  an  empirical  correction  for  the  non-hydrogenic  wave 

functions.  A  description  of  the  method  can  be  found  in  the  original 

29  5  49 

paper  and  in  Griem's  book  ;  a  paper  by  Peach  is  useful  as  an  example  of 

its  application.  Several  comparisons  of  the  predictions  of  this  method  with 

those  obtained  om  more  exact  calculations  using  self-consistent  field 

functions  (see  Section  3.2)  and  with  experimental  results  have  been  per- 

29  4S  50 

formed  by  Burgess  and  Seaton  and  by  Armstrong  and  his  co-workers  ’ 

These  comparisons  show  that  the  method  is  accurate  for  the  photo ionization 

% 

of  an  electron  which  is  alone  in  its  shell  and  in  general  gives  good  agree¬ 
ment  for  light,  more  complex  systems.  Tne  method  of  Burgess  and  Seaton  im¬ 
plies  a  general  restriction  on  photon  energy,  namely  hv-K  <<  z2  Rydbergs, 
this  restriction  is  not  serious  at  the  temperatures  of  interest  in  re  entry 
problems  (<  3  eV)  since  it  occurs  at  photor.  energies  far  from  the  maximum 
of  the  Planck  function.  It  is  nevertheless  worth  recording  that  there 
exists  an  approximation  which  is  valid  for  high  photon  energies  (Ref.  45, 
page  76). 

51  r 

Seaton  gives  a  formula  from  which  -*  may  be  calculated  for  an 

initial  configuration  2p^  (which  corresponds  to  the  most  important  trans^ 
tions  in  nitregeni.  It  is  not  clear  from  the  paper  whether  Seaton's  exprv 
sion  is  obtained  analytically  or  whether  (as  seems  more  likely)  it  is 
empirical  fit  to  other  calculations.  Seaton  suggests  that,  at  least  or 
positive  ions,  the  accuracy  of  the  formula  should  be  within  20%  an  a 
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comparison  due  to  Johnston  et  al.^  with  more  exact  predictions  for  a 
transition  in  Nil  confirms  this  for  low  photon  energies. 

The  final  method  to  be  discussed  is  the  hydrogenic  approximation.  The 
formula  for  the  bound -free  absorption  coefficient  of  a  hydrogenic  gas  is 
(see  e.g.  Kef.  33,  Chap.  5) 


n 


.  !i  e2a2  (l.)3 


3/3 


he 


o  lhw' 


g(v;n,t) 


(5.6) 


n- 


where  hv  is  measured  in  Rydberg  units  and  g(v;n,l)  is  the  Gaunt  factor. 

The  Gaunt  factor  may  be  regarded  as  ?  quantum-mechanical  correction  factor 
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to  the  classical  result,  Defensive  tables  of  Gaunt  factors  are  available 
while  Ref.  53  contains  a  review  of  approximate  analytic  representations; 
except  at  high  photon  energies  the  Gaunt  factor  is  close  to  unity.  The  val¬ 
idity  of  the  hydrogenic  approximation  is  greatest  at  large  values  of  n  and, 
for  given  r.,  for  high  values  of  orbital  angular  momentum.  In  states  of  low 
excitation,  the  core  of  passive  electrons  does  not  fully  shield  the  nucleus, 
and  the  expression  (5.6)  can  be  improved  somewhat  by  the  use  of  zg,  an 
effective  charge,  in  place  of  z.  zg  is  defined  by  z^  s  n2/!^  where 
thr  ionization  energy  1^  is  measured  in  Rydbergs.  When  the  jumping  elec¬ 
tron  is  initially  one  of  a  number  of  equivalent  electrons,  the  hydrogenic 
approximation  breaks  down  completely.  However,  Armstrong  shows  that  the 
p  arentage  splitting  can  be  accounted  for  by  the  fractional  parentage 
coefficient  {F  ( qS * L ' ) } 2  introduced  earlier. 

In  the  nitrogen  calculations  presented  in  Part  3,  bound-free  transitions 
from  the  lower-lying  states  are  represented  by  values  of  cross  section  calcu¬ 
lated  by  Shermian  and  Kulander^0  using  the  Burgess -Sea ton  method.  The  cross 
sections  of  higher  states  are  assumed  to  be  hydrogenic. 
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5,3  Free -Free  Cross  Sections 

Free-free  transitions  involve  the  change  of  energy  of  a  free  electron 
in  the  field  of  a  positive  ion.  When  the  ion  is  hydrogenie  the  following 
expression  holds 


_  16it2e6  2  Ne _  /  kT) 

V’ff  "  3/3  ch  ( 2irm  )3'W/2v*  ’ 


(5.7) 


where  gff(v,kT)  is  the  free-free  Gaunt  factor  which  for  temperatures  of 
interest  to  us,  is  close  to  unity  over  th  important  photon  energy  range 
0.1  eV  to  10  eV;  accurate  values  cf  g_,  ar  tabulated  by  Karzas  and 
umr» 

Until  recently,  the  only  theoretical  treatment  of  free-free  transitions 

55 

was  for  hydrogenic  atoms.  Now,  nowever,  a  paper  by  Peach  is  available 

which  is  an  extension  of  the  Burgess -Seaton  method  (see  Section  5.2)  to 

free-free  transitions.  Peach’s  original  paper  applied  her  method  to  He 

and  compared  the  results  tc  Eq.  (5.7),  differences  up  to  a  factor  of  2  in 

both  directions  occur  with  a  tendency  for  the  hydrogenic  expression  to 

underestimate  the  cross  sections.  Unfortunately  from  our  point  of  view, 

the  results  were  not  averaged  over  electron  velocity  and  it  appears  at 

least  possible  that  this  process  would  decrease  the  discrepancies.  In  the 

2 

past,  it  has  been  common  practice  to  use  an  effective  value  of  z  to 

account  in  a  crude  way  for  non-hydrogenic  effects,  this  has  the  merit  of 

pr  s-^ving  the  analytical  simplicity  of  Eq.  (5.7). 

In  Part  3,  the  free-free  cross  section  is  calculated  according  to  the 

hydrogenic  formula,  Eq.  (5.7).  For  transitions  in  the  field  of  the  first 

2  6l 

ion,  an  effective  value  of  z  based  on  experimental  evidence  A  is  used. 
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APPENDIX  I.  The  Influence  of  a  Density-Dependent  Partition  function  on  the 
Thermal  Equation  of  State  of  Hydrogen 
The  equation  of  state  of  a  gas  is  given  by  the  equation^ 

p  =  MkTV  (I.D 

Then,  using  the  standard  expressions  for  partition  functions  in  a  gas  consist¬ 
ing  of  hydrogen  atoms,  ions  and  electrons  we  get 

p  *  (N„  +  N  +  +  Ng)kT  +  NHVkT  [~  (*n(Qet>H)3T  (1.2) 

H 

which  may  be  written  as  the  sum  of  translational  and  electronic  contributions 
to  the  pressure, 

P  ’  Ptr  +  Pet-  <I-3) 

We  see  from  Eq.  (1.2)  that  the  electronic  partition  function  can  affect  the 
value  of  the  pressure  in  two  ways.  In  the  first  p.Taco  it  determines,  through 
the  Saha  equation,  the  number  of  particles  present.  The  particle  densities 
can  be  accurately  calculated  in  our  range  cf  conditions  ii'  the  truncation 
point  of  the  partition  function  is  found  from  Debye-Htlckel  theory  (see  Chap.  2). 
This,  however,  leads  to  the  second  way  in  which  the  partition  function  influ¬ 
ences  the  pressure.  The  use  of  Debye-HUckel  theory  leads  to  a  density- 
dependent  partition  function  and,  according  to  Eq.  (1.2),  this  gives  rise  to 
a  cnange  in  the  thermal  equation  of  state  in  the  form  of  an  additional  term 
which  we  call  p  in  Eq.  (1.3).  It  is  this  additional  term  which  v»  discuss 
now. 

The  terms  of  the  hydrogenic  partition  function  are  conveniently  ordered 
by  the  principal  quantum  number  n  and  we  may  write 
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‘  max  -E  /kT 

°et  "I  *n  *  " 


(1.4) 


where  g  s  2n2  and  E  =  xu(l"l/n2)*  Following  the  procedure  of  Appendix  II 
n  n  h  * 

we  may  separate  Qel  into  two  parts, 


Qei  s  (Qet'low  +  *Qei*high 


(1.5) 


where 


n  -1  -E  /kT 

(Q  i  I  2n2  e  n 

et  low  L. 


(1.6) 


(n  ,  s  -E  /kT 

^et  high  “  l  *  2n2  e  n 


-E  /kT 

2n2  e  n  dn. 


(1.7) 


We  now  choose  n  such  that  the  error  in  replacing  E  by  x  u  is  small  and 

n  h 

we  get 


-X  /kT 

“Vhigh  5  2/3  e  "  «V»>’  -  ("V]. 


(1.8) 


Now  since  ^e£^ow  independent  of  V,  we  get  from  Eqs.  (1.2),  (1.5)  and 

(1.8). 


,y2  -XH/kT 

n  a  N  VkT  — —  e  "  (  _) 

Pel  H  Q  £  3V  T’ 


(1.9) 


Where  we  have  replaced  n^  by  v  to  emphasize  that  we  are  treating  it  as 
a  continuous  variable,  v  is  found  from  the  Debye-HUckel  theory  discussed  in 


60 


Section  2.2,  Eqs.  (2.6)  and  (2.7),  together  with  the  hydrogenic  expression  for 
for  energy,  Cq.  (2.S),  those  combine  to  give 


XH  kT  1/,( 

V  .  {JL— EL_)  . 

e4  8we2H 


(1. 10) 


Writing  a  for  the  degree  of  ionization  so  that  a  = 

we  get 


v(i vy  *  H 


.  .  (.S.  — SI _ X.) 

e4  8xe22oNL  Vo 


(1. 11) 


so  that 


.  1  1 

“  4  V  * 


(1.12) 


Substituting  Eq,  (1.12)  into  Eq.  (1.9)  we  find 


1  *  -XH/kT 

*  2  V  v  •  1/Qef 


(1.13) 


Since  the  translational  pressure  due  to  the  hydrogen  atoms  is  given  by 

(pu).  1-  N  kT,  this  may  be  written  as 

n  tp  n 


1  ,  "XH/kT 

P.(  a*2  Vtr  v  *  1/Q, 


(1.14) 


Now,  identifying  n  with  v  and  comparing  Eq.  (1.8)  we  see  that  the 

ITI3X 

1  ,  "XH/kT 

quantity  ^  v3  e  is  related  to  the  high  state  contribution  to  the 

partition  function  and  we  make  the  definition 


_  i  _XH^k'^  3  * 

5  5  2 v  •  *  r  C(Vhigh  * 0 1 


(1.15) 
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where 


*  *  ,  -XH/kT 

Q  5  2/3(n  )3e  H 


(1.16) 


In  addition  we  have 


ai  ^  ‘’it  « 

tr  _ H _ _  1  -  a 

ptr  ",  +  M...  +  NH  1  +  0 


(1.17) 


so  that  Eq.  (1.14)  finally  becomes 


_®i  -  _1_  1  ~  « 

’’tr'O.l1’0  ‘ 


(1.18) 


For  v  to  be  approximately  continuous  as  required  by  the  derivation 
of  Eq.  (1.18),  its  value  has  to  be  relatively  large  which  implies  low  den¬ 
sities  which  in  turn  implies  a  high  degree  of  ionization.  Equation  (1.18) 
shows  that  under  these  circumstances  Pe£/Ptr  will  b®  small.  Numerical 
checks  confirm  that  for  a  >  0.85,  Pe£/Ptr  <  °*°5.  The  physical  interpre¬ 
tation  is  that  there  are  so  few  atoms  present  that  the  change  in  their 
partition  function  is  irrelevant  to  the  total  pressure. 

Where  the  ionization  is  low,  v  tends  to  be  low  (because  the  density 
is  usually  high)  and  the  integral  approximation  to  ^e£^igh  *8  *688 
appropriate.  However,  the  error  is  probably  no  worse  than  that  caused  by 
treating  the  energy  levels  as  unperturbed  (which  is  the  only  treatment 
available  at  present).  We  therefore  retain  the  integral  approximation  for 
small  values  of  v.  Thus,  Eq.  (1.18)  is  assumed  to  hold  but  now  we  need 
to  evaluate  $/Qel  more  carefully,  replacing  Q  by  its  definition  in 
Eq.  (1,15)  ve  obtain 


•'et  1  e 

*  2  Q 


» .V“ 


1-0 
1  ♦  o 


(1.19) 
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Numerical  evaluation 
5%  for  densities  up  to  1 ) 


of  Eq.  (1.19)  shows  that  p  ,/p^  is  less  than 

e*  tr 

snd  temperatures  up  to  3eV. 
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APPENDIX  II.  The  Partition  Functions  of  Nitrogen  and  its  First  Ion 

This  appendix  is  primarily  concerned  with  giving  details  of  the  energy 

levels  and  statistical  weights  from  which  the  partition  functions  shown  on 

Figs.  2  and  3  were  calculated.  Most  of  the  low  and  moderate  energy  levels 

57 

were  obtained  from  the  compilation  of  Gilmore.  This  compilation  conven¬ 
iently  combines  closely  spaced  energy  levels  of  each  electronic  configu¬ 
ration.  The  basic  data  for  Gilmore's  tabulation  consisted  of  exnerimental 
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results  listed  by  Moore  supplemented  with  Gilmore's  own  estimated  values. 

Gilmore's  work  covers  states  with  principal  quantum  numbers  up  to  8.  Above 

58 

this  level,  Moore's  values  were  taken  and  combined  where  available  (the 
method  of  combination  is  given  in  a  later  paragraph).  Otherwise,  hydrogenic 
estimates  asymptotic  to  the  true  ionization  thresholds  were  made.  For  very 
high,  closely  spaced  energy  levels  integral  approximations  to  the  summed 
terms  were  used. 

The  electronic  partition  functions  were  calculated  from  expressions  of 
the  type 


Q«1  V 

The  meaning  of  the  terms  is  as  follows 

i(D 

S1  2  fX  «[])  **p("  Ei1>/kT) 

1  i*l  1  1 

and  is  a  sum  over  all  energy  levels  with  n  3.  These  levels  are  below 
the  lowest  value  to  which  the  ionization  energy  is  depressed  in  our  range 
of  conditions.  This  means  that  the  associated  electronic  energy  configure 
tion  of  each  level  need  not  bo  distinguished  and  occasional  further  combi¬ 
nation  of  Gilmore's  values  has  been  possible. 
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and  Sg  are  defined  in  a  similar  manner. 


(2) 

"max 


c,  =  I  g<2)  exp<-  E^VkT) 

-  i=l  1 


,<2) 


and 


.(3) 

i 

max 


s,  5  l  k{3>  «tp(-  E^'/kT), 
i*l 


.(3) 


These  sums  cover  quantum  numbers  4  <  n  <  10.  Zr  :h  corresponds  to  a  differ¬ 
ent  electronic  core  configuration.  This  distinction  is  necessary  because 
for  a  given  depression  in  ionization  potential.  Ax,  the  cut-off  energy 

(x  s  X~6x)  depends  upon  the  ionization  energy  and  this  takes  different 
«rr 

values  according  to  the  core  configuration.  A  word  about  the  values  listed 
in  Tables  II. 1  and  II. 2  is  necessary.  Some  of  the  energy  levels  of  and 

Sg  could  have  been  combined  without  significant  loss  of  accuracy  but  were 
individually  specified  for  aesthetic  reasons.  The  levels  in  question  tend 
to  have  large  valuer  o*  statistical  weight  and,  when  the  series  cut-off  is  in 
this  i'egion,  the  function  suffers  large  discontinuities  as  a  function 

of  density  (because  we  neglect  perturbation  of  the  energy  levels).  These  dis¬ 
continuities  lead  to  nonsmooth  functions  for  the  particle  densities.  This 
lack  of  smoothness  can  be  avoided  by  including  a  sufficient  number  of  terms 
in  the  sums  and  S^. 

The  quantities  1^  and  I3  are  integrals  defined  as  follows 


max 


I„  2  g<2>  «*P<  rl:  *T>  2n2or 


’cor* 


11 


*  2,3  *111.  EiJ)/kT)(nL,  -  1331> 


and 
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:3  2  Score  exp(-  Ei3>/kT)  {  'WX  *>’<>» 


■■  5/3  <£r.  «P<‘  Ei3)/kTXn3  -  1331). 


(2)  (3) 

where  and  g _  are  the  statistical  weights  of  the  corresponding 

core  core 

core  configurations.  and  1^  are  therefore  integrals  over  high  states 

whose  energy  levels  are  approximated  by  the  threshold  values. 

Values  of  the  energies  and  statistical  weights  are  given  for  the  atom 
in  Table  II. 1  and  for  the  first  ion  in  Table  II. 2.  It  retrains  to  demon¬ 
strate  the  method  of  combining  closely  spaced  energy  levels. 

For  q  closely  spaced  energy  levels  we  require  a  mean  energy  level  £ 
such  that 


q  -Ej/kT  -E/kT  q 

l  «  *  •  Jr* 

r-i  *  i-s  i 


(ii. l) 


If  we  now  write  *  g-fAEi  and  substitute  in  Eq.  (II. 1)  we  obtain 


q  -AEj/kT  q 

l  *i  •  *  I  Ri 

i*l  1  i*l  1 


(II. 2) 


Sj*c  if  AE^'kT  «  lt  Eq.  (TL2)  becomes 


j.  *i“i  * 0 


(II. 3) 


which  yields 


1  *ici 


i»l  * 
9 

l  R* 
»  _*  • 


(II.4) 
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Equation  (II. 4)  can  therefore  be  used  to  combine  any  number  of  energy  levels 
provided  only  that  in  each  case  |E.-E|  «  kT. 


TABLE  II.l 


Energy  Levels  and  Statistical  Weights  of  N 
Energy  Values  in  Electron  Volts 


Core: 

2s22p2(3P) 

2s22p2(‘D) 

i 

,  (1> 
*i 

'itl) 

(2) 

*i 

E.(}> 

1 

„  <3>  P  O) 

*i  Ei 

1 

4 

0 

54 

13.318 

160  15.448 

2 

10 

2.384 

216 

13.686 

250  15.903 

3 

6 

3.575 

450 

13.957 

360  16.069 

4 

18 

10.450 

648 

14.158 

490  16.169 

5 

12 

10.926 

882 

14.270 

• 

640  16.234 

6 

54 

11.875 

xi52 

14.335 

810  16.279 

7 

10 

12,356 

14  58 

14.380 

1000  16.311 

8 

108 

12.980 

1800 

14.412 

9 

30 

13.758 

10 

66 

15.093 

Core 

statistical  weights 

(2) 

*core 

*  9 

(3)  e 

g  *5 

core 

Ionization  energies 

*L7) 

*  14.5481 

Ei3)  *  16.447 

*  16.447 
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TABLE  I I. 2 


Energy  Levels  and  Statistical  Weights  of  N' 


Energy  values  in  electron  volts 


Core: 

2s22p(2P) 

2s2p2(4P) 

• 

1 

(1) 

«i 

E.(1) 

l 

,  <2> 
gi 

e.<2> 

1 

P  (3) 

gi 

e.<3) 

1 

1 

9 

0 

192 

25.860 

384 

32.990 

2 

5 

1.899 

300 

27.337 

600 

34.537 

•> 

w 

1 

4.053 

432 

28.081 

864 

35.201 

4 

5 

5.800 

,8C 

28.500 

1176 

35.601 

5 

15 

11.436 

768 

28.7S0 

1536 

35.861 

6 

9 

13.541 

972 

28.940 

1944 

36.041 

7 

5 

17.876 

1200 

29.068 

2400 

36.169 

8 

12 

18.480 

9 

3 

19.232 

10 

39 

20.933 

11 

60 

23.270 

12 

221 

29.045 

Core 

statistical  weights 

(2) 

£  = 
^core 

6 

(3) 

e 

^core 

=  12 

Ionization  energies 

e(2)  * 

• 

29.612 

e(3) 

=  36.713 
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APPENDIX  III.  The  Approximations  Involved  in  the  Use  of  a  Density- 
Independent  Partition  Function 

In  radiative  transfer  calculations  we  require  the  species  composition 
in  order  to  find  the  occupation  numbers  of  the  initial  levels  of  the  various 
radiative  transitions  and  we  also  require  the  electron  and  positive  charge 
densities  to  find  the  widths  of  the  stark  broadened  lines.  It  turns  out 
that  these  quantities  are  not  very  sensitive  to  the  cut-off  of  the  partition 
functions  and  can  be  calculated  reasonably  accurately  without  allowing  for 
the  high  terms  of  the  partition  functions.  This  appendix  discusses  the 
effect  on  occupation  numbers  and  charged  particle  densities  of  truncating 
the  partition  function  after  a  small,  fixed  number  of  terms. 

The  possibility  of  such  an  approximation  springs  from  the  following 
observations.  In  the  first  place,  at  low  temperatures  the  factor  1/kT 
in  the  exponents  of  the  terms  of  the  partition  functions  means  that  the 
high  terms  do  not  contribute  significantly  and  the  cut-off  is  unimportant 
(except  for  its  existence).  (From  this  fact  it  follows  that  the  problem  is 
restricted  to  relatively  high  temperatures  where  the  gas  is  completely  dis¬ 
sociated  and  we  ignore  the  molecular  species  in  the  remainder  of  the 
discussion.)  Secondly,  excited  states  of  multiply  ionized  species  have 
large  energies  with  respect  to  the  ground  state  and  the  corresponding  con¬ 
tributions  to  the  partition  functions  are  small,  the  multiple  ions  thus  tend 
to  be  dominated  by  a  few  low-lying  states  and  their  partition  functions  are 
consequently  density  independent.  Checks  performed  a  posteriori  show  that 
we  do  not  need  to  consider  the  problem  of  the  exact  truncation  point  for 
any  multiply  Ionized  species  in  our  range  of  conditions.  Furthermore, 

although  (Q  t)  is  treated  as  density  dependent  in  the  accurate  calcula- 
N 

tions,  it  does  not  vary  more  than  a  factor  of  3.5  from  a  sum  over  its  first 
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four  terms  (see  also  Fig.  3).  Thus  It  can  be  seen  that  only  (0  „),.  is 

et  N 

strongly  dependent  on  density. 

It  follows  that  when  few  atoms  are  present,  the  electron  and  positively 
charged  particle  densities  may  be  calculated  reasonably  accurately  no  mat- 
ter  what  form  of  the  atomic  partition  function  is  used.  On  the  other  hand, 
a  significant  number  of  atoms  (i.e.  low  ionization)  will  only  be  present  if 
the  temperature  is  low  or  the  density  is  high.  If  the  temperature  is  low, 
the  terms  emitted  from  the  approximate  partition  function  are  negligible, 
as  explained  previousl>  If  the  density  is  high,  there  will  be  a  large  de¬ 
pression  in  the  ionization  limit  and  thus  the  density-dependent  cut-off 
point  will  be  relatively  close  to  the  low-lying,  density-independent  value 
which  will  now  be  a  better  approximation.  In  summary,  if  we  approximate 
the  atomic  partition  function  by  its  first  few  terms,  then  under  conditions 
of  low  ionization,  where  the  value  of  the  atomic  partition  function  is  im¬ 
portant,  the  approximation  is  reasonabxy  good  while  for  conditions  of  high 
ionization,  where  the  approximation  is  very  poor,  the  charged  particle 
density  is  independent  of  its  value. 

At  conditions  where  the  ionization  is  low,  the  approximate  partition 
function  is  relatively  accurate,  as  discussed  above,  we  can  therefore  use 
it  to  find  the  occupation  numbers  of  the  atomic  states  with  reasonable  accu¬ 
racy.  There  remains,  however,  the  problem  of  their  calculation  for  higher 
ionization  (where  the  atomic  radiative  contribution  is  not  necessarily 
negligible).  Under  these  conditions,  a  further  observation  is  of  help. 
Although  the  region  of  high  ionization  is  the  region  of  largest  error  in 
the  density- independent  partition  function  (see  Fig.  2),  it  is  also  a  region 
in  which  the  occupation  numbers  become  independent  of  the  partition  function. 

This  may  be  seen  by  noting  that  N  and  N  will  be  independent  of  NVI 

e  Nt  N 

and  using  Eq.  (2.1)  to  eliminate  NN  from  the  appropriate  form  of 
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Eq.  (2.2),  to  obtain 


N  N 

(T~T~  nH  (kT>  0/po  *i  expC"  (14*548-AX1-E.)/kT3  (III.l) 

ei  Nt 


59 

A  similar  result  has  been  independently  obtained  by  Hochstim.  The  value 
of  is  exactly  linear  in  under  these  conditions  and  the  influ¬ 

ence  of  inaccuracies  in  (Q  tf)M  can  be  interpreted  by  saying  that  although 

f5  it  n 

we  may  not  know  accurately  the  total  number  of  atoms  present  —  because  we 
do  not  know  whether  to  count  highly  excited  states  as  bound  or  free  —  once 
having  decided  to  treat  a  particular  state  as  bound  we  can  calculate  the 
associated  occupation  numbers  with  reasonable  accuracy. 

The  solution  of  the  system  of  equations  consisting  of  the  first  four 
ionization  equations,  Eq.  (IV. 12)  to  (IV. 15)  plus  the  two  conservation  equa¬ 
tions,  Eqs.  (IV. 16)  and  (IV. 17)  was  obtained  by  iterative  calculations  on 
an  IBM  360  computor  using  both  exact  and  approximate  expressions  for  (Qei)N 

and  (Q  )  ,  The  approximate  expressions  consist  of  a  fixed  number  of  the 

e  N 

terms  listed  in  Appendix  II  while  the  "exact’  calculations  are  described  in 
Section  2.2. 

The  composition  predicted  by  the  "exact"  theory  is  shown  on  Figs.  4-7. 
Figures  15,  16  and  17  display  some  of  the  errors  incurred  by  using  vacuum 
ionization  energies  and  the  following  truncated  partition  functions  for  N 
and  N+ 


(Qe£)N  =4+10  exp(-  2.384/kT)  ♦  6  exp(-  3.575)  (III. 2) 

<Q  .)  *=9+5  exp(-  1.899/kT )  ♦  exp(-  4.053/kT)  ♦  5  exp(-  5.6uO/kT) . (HI.3) 
et  N+ 

Figure  15  shows  the  errors  in  electron  number  density  which  for  mass  densities 
less  than  1  atmosphere  are  within  15%.  Because  the  overall  charge  must  be 
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zero,  the  total  density  of  positive  charges  must  equal  the  electron  density 
and  does  not  need  to  be  considered  further.  The  density  of  atoms  suffers 
the  worst  errors,  as  is  to  be  expected,  and  Fig.  16  shows  that  these  can  be 
two  or  three  orders  of  magnitude;  N*  is  only  seriously  in  error  when  mul¬ 
tiple  ionization  dominates  (Fig.  17). 

The  error  in  the  occupation  numbers  for  various  approximations  is 
shown  in  Figs.  18-20.  Figure  18  is  based  on  the  3-term  partition  function 
of  Eq.  (III. 2)  and  vacuum  ionization  potentials.  The  maximum  error  occurs 
when  the  atoms  are  60%-70%  ionized  which  means  that  it  moves  towards  higher 
densities  as  the  temperature  rises,  and  it  grows  from  2%  at  1.0  eV  to  over 
60%  at  3  eV.  For  densities  less  than  1  atmosphere,  the  error  is  less  than 
40%.  Figure  19  shows  the  errors  caused  entirely  by  the  partition  function 
(the  calculations  plotted  on  Fig.  19  used  correctly  depressed  ionization 
energies  but  the  simple  3-term  partition  function  of  Eq.  (III. 2)). 

The  fact  that  the  errors  are  worst  at  high  densities  and  that  the 
partition  function  and  ionization  errors  are  both  in  the  same  direction  sug¬ 
gests  that  one  may  be  able  to  reduce  the  errors  and  perhaps  compensate  for 
the  heqvily  reduced  ionization  energy  at  high  densities  by  truncating  the 
partition  function  slightly  above  the  exact  high-density  cut-off.  Figurt  20 
shows  the  results  of  such  a  procedure  where  the  exit -off  point  of  the  atomic 
partition  function  was  chosen  to  be  the  true  cut-off  for  densities  slightly 
less  than  atmospheric.  This  gives  rise  to  a  13-term  partition  function  of 
the  form 


l  «!1>  «*p(-  E(1>/kT)  ♦  l  g[2)  exp(-  E(2)/kT) 
i*l  1  i*l 

♦  g^  exp(-  E^VkT). 


(III. 4) 
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The  partition  function  of  N+  was  left  in  its  approximate,  4-term  form. 
The  maximum  error  can  be  seen  from  Fig.  20  to  have  been  reduced  to  20%. 
However,  the  choice  of  cut-off  turns  out  to  be  critical  since  use  of  the 
true  cut-off  corresponding  to  10”^  atmospheres  increases  the  high  density 
error  to  -50%.  It  also  turns  out  that  severe  errors  can  occur  in  the 
occupation  numbers  of  N+  when  the  atomic  partition  function  is 
overestimated. 

In  conclusion,  an  approximate  partition  function  based  on  a  small 
number  of  low-lying  energy  levels  gives  occupation  numbers  and  charged 
particle  densities  to  an  accuracy  of  about  50%  provided  the  densities  are 
less  than  atmospheric.  In  many  cases,  this  will  be  an  acceptable  approxi¬ 
mation.  However,  attempting  to  reduce  the  inaccuracies  by  raising  the 
cut-off  point  of  the  partition  function  introduces  uncontrolled  errors  and 
does  not  seem  to  be  a  fruitful  approach. 
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APPENDIX  IV.  An  Iterative  Procedure  for  the  Solution  of  the  Mass  Action  and 
Conservation  equations 

At  high  temperatures,  a  diatomic  gas  A0  will  undergo  dissociation, 


A2  -  2A, 


(IV. 1) 


ionization  of  the  molecules. 


>*-l  T  .z 


a2  ♦  •  , 


(IV. 2) 


(where  z  denotes  the  number  of  positive  charges  on  the  molecular  ion)  and 
ionization  of  the  atoms. 


z-1  -*■  z 
A  ♦  A  ♦  e 


(IV. 3) 


The  corresponding  equations  of  mass  action  can  be  obtained  from  Eq.  (2.2)  by 
identifying  the  symbols,  inserting  the  standard  approximation 

(Q)g  =  (Qtr^s^int an<i  evaluatin*  ^tr^S*  The  resultinR  equations  are, 
for  dissociation: 


‘V*.,’  -DAT 


(IV. 4) 


where  D  is  the  dissociation  energy,  is  the  mass  of  an  atom  of  A 

and  i*  the  electronic  partition  function  of  A;  for  molecular 


ionization: 


N  N  (0.  > 

A*  *  int  A*  2w  a  kT  3/2  XA*'X 

2  -  2  *  2  *  ( — 2-)  «xp(-  -Ir-) 


n  ,  *  To TT  . 

.z-1  int  .z-1 


h2 


(IV. 5) 


?«  m  3/2 

where  subscript  e  denotes  electrons  and  the  value  of  ( - )  is 

3.03  -  1021  eV~3/2  cm"3; 


and  for  atomic  ionization: 
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N 


N 


tj. .  _'V*L .  2 ,  f2’ 

r1  (0et)Az-l  1,2 


3/2 

)  exp(-x  ,/kT).  (IV. 6) 

Az-i 


Note  that  it  is  possible  to  reorganize  the  sets  of  equations  (IV.4),  (IV. 5) 
and  (IV. 6)  so  that  molecular  ions  are  not  regarded  as  being  produced  by 
ionization  but  instead  treated  as  dissociating  according  to  an  equation 
similar  to  Eq.  (IV.u);  this  is  exactly  equivalent  to  the  present  treatment 
and  the  further  inclusion  of  dissociation  relations  for  molecular  ions 
would  introduce  redundant  equations. 

The  conservation  equations  for  mass  and  charge  are 


2 


2z 


I  " ,  ♦  l  " 


z*0  A 


z*0  A 


(IV. 7) 


and 


a  a 

l  *N  ♦  l  zN  •  N  (IV. 8) 

*«1  A*  z«l  A* 


where  N  is  constant  and  is  the  number  of  atonic  nuclei  per  unit  volume 
o 

and  is  the  atomic  number  of  A.  If  a  method  of  evaluating  the  parti* 
tion  functions  is  adopted,  the  set  of  equations  (IV.4)  to  (IV. 8)  can  be 
solved  by  iteration. 

In  the  case  of  nitrogen,  the  only  molecular  ionisation  process  of 
importance  is  the  reaction 

Nj  *  N*  ♦  e"  (IV. 9) 

and  further  it  can  easily  be  checked  that  for  kT  «  3  eV  and  o/o  >  10*^. 

o  * 
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the  fifth  and  higher  ions  never  contribute  more  than  1%  of  the  total 
particle  density  and  can  be  ignored.  We  therefore  may  truncate  the  summa¬ 
tions  in  Eqs.  (IV. 7)  and  (IV. 6)  after  z  «  1  for  molecules  and  z  *  4  for 
atoms |  in  addition  we  have  to  consider  one  equation  of  type  (IV.S)  and  four 
of  type  (IV. 6). 

Taking  account  of  the  simplifications  discussed  above,  the  explicit 
form  of  equations  (IV. 4)  to  (IV. 8)  for  nitrogen  may  be  written  as 


V  ".  ♦ 

1/9  .1 

- -  K.  *  112.4  ...  ■■£  (k?r"  (p/p  )  A  exp(-  15.6/kT)  (IV.10) 

W„  1  'QintX  ° 


U2  .  3/2 

—  «  ■  0.81  *  108  (kT)3'2 


int  H. 


(p/p  )  exp(-  9.56/kT)( IV. 11) 


N  N  (Q  ,) 

u*  3/2  ,i 

-2 - -  K  «  112.4  ——JL.  (kT )*'*  (p/p  )  A  exp[-( 14. 548-Ax. )/kT] 

M..  '"et'N 


(IV. 12) 


N  N  (Q) 

3/2  -1 

\  •  112.4  -  (kT)  "  (p/po)  A  exp(-( 29.605-AXj )/kT ] 

Mm^  (IV.  13) 


- —  ■  Kj  •  112.4  p -  (kTr"  (p/po)  A  expC-(47.463-Ax3)/kT) 

*  *l  N**  (IV. 14) 

n 


N  A  H  (Q  ,)  _ 

„♦♦♦♦  e  el  ♦♦♦♦  3/2  , 

-2 - -  K  *  112.4  — - J: -  (kT)  '  (p/p  )  exp[<T;.45-Ax4)/kT] 

W...  %*„♦♦♦ 


(TV. 15) 
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2$  ♦  2NV,  ♦  N„  ♦  N  ♦  M 


N 


N 


.♦t 


♦  N 


♦  N 


(IV. 16) 


N  ♦  N  ♦  2N  t  3N  x  4  4H  * 

N*  „♦  H++  H+++  if*** 


e 


(IV. 17) 


where  the  tilde  denotes  a  particle  density  non-dimensional ized  by  KQt  the 

total  number  of  nuclei  present  and  Axt  is  the  ionisation  depression  of 

Eq.  (2.7),  The  reductions  in  ionization  energies,  (0<£)N  and  (0el>  ^ 

N 

are  all  density  dependent  and  depend  for  their  evaluation  on  a  previous 
solution  of  the  mass  action  equations,  as  has  been  described  earlier;  in 
this  appendix,  however  we  regard  them  as  given  and  describe  the  iterative 
solution  of  Eqs.  (IV. 10)  to  (IV. 17). 

The  solution  of  this  set  of  equations  is  made  considerably  easier  by 
the  fact  that  at  most  two,  and  frequently  only  one,  of  the  species  are  much 
more  abundant  than  the  others,  this  can  be  seen  as  follows.  For  i  ^  2  we 
have  <<  K^,  consider  first  the  case  where  few  uncharged  particles  are 

present  so  that  1  <  N#  <  4  then,  if  N*  is  a  x- times  ionised  particle 


s 

and  the  species  N  °  and 
will  dominate  since 


ii  .  •  H  K 

H*  **2 


*o  *  for  which  K  _ 
It  s  *2 

o 


is  nearest  to  unity 


and 


«  it 

H 


s  *1 
o 


for  s  >  s  ♦  2 
—  o 


<< 


for  s  <  z  -  1 
—  o 
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by  virtu*  of  th«  values  of  K  ,  In  th*  case  where  there  are  a  significant 

•  \  « 

number  of  uncharged  particles  present,  N  *  N  and  it  again  turns  out 

*  N* 

that  the  species  for  which  *8  near  unity  will  dominate.  N  +  can 

never  be  Important  however,  since  although  can  be  unity,  implying 
M  «  av  the  same  time  K.  >>  1  so  that  H  <<  (N  )2.  The 

Hj  *2  2  "j  " 

criteria  for  picking  the  dominant  equation  are  not  critical  since  although 

it  can  happen  that  «  1  and  >>  1  such  that  K  *  1  (i.e. 

H  .  •  M  ,)  at  th*  same  time  N  ,  >>  M  .  *  M  .  and  either  equa- 

M*’1  N*'3  K*“2  K*"1  N*“3 

•  e»2 

tioo  will  give  the  dominant  species  N  .  A  detailed  consideration  of  the 


equations  leads  to  the  following  criteria  for  the  dominant  equation 


Kj  <fr3  «  1, 


Eq.  *IV.ll)  is  dominant 


Kj  d?  >  1  and  ^  1,  Eq.  (IV.  12)  is  dominant; 


>  1  and  4,  Eq.  (IV. 13)  is  dominant; 


>  *  and  *5K6  <_  9,  Eq.  (IV.  14)  is  dominant 


K,K  >  Eq.  (IV.  15)  is  dominant. 


Since  w*  have  one  dominant  equation  which  fer  given  conditions  we  can 
identify,  we  can  solve  the  set  of  equations  by  an  iterative  process  organ¬ 
ised  as  follows.  Suppose,  first,  that  the  dominant  equation  is  one  of 
Eqe,  (IV. 12)  to  (IV. 15).  It  therefor*  has  th*  form 
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* 

N 

N 


z+1 


N 

e 


K 


(IV. 18) 


and  Eqs.  (IV. 16)  and  (IV. 17)  may  be  written 


ft  ♦ 
z 


N 


*♦1 


(IV. 19) 


and 

zN  ♦  (ztl)ft.  *  ft  -  ft  .  ( IV. 20) 

z  ztl  e  c2 

reapectivaly  where  z  >  0  and  ft^,  ^c2  ar*  the  terms  in  Eqs.  (IV. 16)  and 
(IV. 17)  not  individually  exhibited  in  Eqs.  (IV. 19)  and  (IV. 20)  and  are 
regarded  as  small  corrections.  Eliminating  ft  ^  from  Eq.  (IV. 20)  the 
equations  can  be  written  in  the  form 


ft  «  A  -  ft  (IV. 21) 

z+1  z 

ft  *B-ft  ■  B  -  A  ♦  ft  .  (IV. 22) 

e  z  ztl 


where 

A  s  l  -  ft  (IV. 23) 

cl 

B  I  (z+l)A  ♦  ft  _  (IV, 24) 

C2 

If  we  regard  A  and  B  as  known,  Eqs.  (IV. 18),  (IV. 21)  and  (IV. 22) 

may  be  easily  solved.  Tree  the  resulting  values  of  N  and  N^,  the 

particle  densities  of  the  remaining  species  may  be  calculated  by  Eq.  (IV. 10) 

and  the  four  subsidiary  equations  out  of  (IV. 11)  to  (IV, 15),  corrected  values 

of  N  ,  and  H  _  can  then  be  calculated  and  a  higher  approximation  to  the 
cl  c? 
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to  the  dominant  species  obtained.  Checks  show  that  if  the  iteration  is 
continued  until  N  is  within  0.5%  then  all  species  are  within  0,05%.  If, 
now,  the  dominant  equation  is  Eq«  (IV. 11),  a  similar  procedure  can  be  used 
although  the  equations  are  slightly  different  in  detail.  We  have 


(»■>* 


N, 


k2. 


N, 


(IV. 25) 


and 


NN  ■  A  ' 


(IV. 26) 


where 


A  »  1  -  N  ,  *  1  -  (2H  ♦  N  ♦  N  -t  N  ♦  N  )  (IV. 27) 

el  H*  H+  Ht++ 


Tnis  iterative  procedure  is  rapidly  convergent  and  has  been  programmed  as 


a  FORTRAN  subroutine 
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FIG  10  SKETCH  SHOMflNG  A  SIMPLE  ARRANGEMENT 
OF  ENERGY  LEVELS  ANO  ILLUSTRATING 
SOME  TYPICAL  BOUND  -  BOUNO  TRANSITIONS 
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OF  ENERGY  LEVELS  AND  ILLUSTRATING 
SOME  TYPICAL  BOUND -FREE  AND  FREE  - 
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FIG.  4  PARTICLE  DENSITIES  IN  A  NITROGEN  PLASMA 
AS  A  FUNCTION  OF  DENS'TY  FOR 

*T  2  0.5  eV  (Ta  5,803  °K  ) 

N0*5.374  x  io  9  xp/pQcm  5, /O0-i. 250*10""  3gm  cm 


FIG. 5  PARTICLE  DENSITIES  IN  A  NITROGEN  PLASMA 
AS  A  FUNCTION  OF  DENSITY  FOR 
kT -I  OeV  (Ts  llt605  PK  ) 

V5  :574x,°,9x/)/iP0cm'“5,p0=|.25rsviO'3gm  cm'3 
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FIG.  6  PARTICLE  DENSITIES  IN  A  NITROGEN  PLASMA 
AS  A  FUNCTION  OF  DENSITY  FOR 
kT  =  2.0 eV  (Ts  23,211  °K) 

N0  =  5. 374X10  xp/po cm  i/>o  =  l.250x|0  gm  cm 
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FIG.  7  PARTICLE  DENSITIES  IN  A  NITROGEN  PLASMA 
AS  A  FUNCTION  OF  DENSITY  FOR 

kT  =  3.0  eV  (Ts  34,8I6°K) 

N0=5.374x|Ol9xn/o0cm“3l/?0=l.250XIO’3gm  cm“ 


FIG.  80  ABSORPTION  COEFFICIENT  FOR  A  LINE 
WITH  A  DISPERSION  PROFILE. 


{i-exp{-[ni- 


AS  A  FUNCTION  OF  FREQUENCY  FOR  A 
TYPICAL  SELF-  ABSORBED  LINE. 


FIG.  8C  EQUIVALENT  WIDTH  AND  SPECIFIC  INTENSITY 
AS  A  FUNCTION  OF  FREQUENCY  FOR  A 
TYPICAL  NONABSORBED  LINE. 


ALL  LINES  ^ 
EVENLY  SPACED 


FIG.  10  A  SKETCH  OF  POSSIBLE  CURVES  OF  GROWTH  OF  A  GROUP  OF 
LINES.  (THERE  IS  NO  SIGNIFICANCE  IN  THE  RELATIVE  POSITIONS 
OF  THE  TRANSITION  REGIONS  ON  THE  Tm  SCALE) 
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FIG.  II  THE  CURVE  OF  GROWTH  OF  A  PURE  DOPPLER  LINE 
PROFILE  AND  ITS  APPROXIMATE  REPRESENTATIONS 
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FIG  12  NORMALISED  LINE  WIDTHS  OF  ELECTRON  IMPACT  BROADENED 
PROFILES  IN  THE  FIRST  SPECTRUM  OF  NITROGEN.  (VALUES 
TAKEN  FROM  GRIEM  W) 


FIG  13  NORMALISED  LINE  WIDTHS  OF  ELECTRON  IMPACT  BROADENED 
PROFILES  IN  THE  SECOND  SPECTRUM  OF  NITROGEN  (VALUES) 
TAKEN  FROM  GRlEM0j) 
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~IG.  14  CURVES  OF  GROWTH  OF  COMBINED  DOPPLER  AND 

DISPERSION  LINE  PROFILES  BASED  UPON  THE  CORRESPONDING 
PURE  DISPERSION  PROFILE  PARAMETERS. 
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FIG.  15  ERRORS  IN  ELECTRON  DENSITY  DUE  TO  USING  VACUUM 

IONISATION  ENERGIES.  A  3-TERM  PARTITION  FUNCTION 
FOR  N  ANO  A  4-TERM  PARTITION  FUNCTION  FOR  N  + 
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FIGS  16  ond  17  ERRORS  IN  NUMBER  DENSITIES  OF  N  AND 

N*  DUE  TO  USING  VACUUM  IONISATION 
ENERGIES,  A  3-TERM  PARTITION 
FUNCTION  FOR  N  AND  A  4-TERM 
PARTITION  FUNCTION  FOR  N* 
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FIG  18  ERRORS  IN  THE  OCCUPATION  NUMBER  OF  THE  iTH 
ENERGY  LEVEL  OF  N  DUE  TO  USING  VACUUM 
IONISATION  ENERGIES,  A  3-TERM  PARTITION 
FUNCTION  FOR  N  AND  A  4-TERM  PARTITION 
FUNCTION  FOR  N  * 
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FIG.  IS  ERRORS  IN  THE  OCCUPATION  NUMBER  OF  THE  iTH  ENERGY 
LEVEL  OF  N  DUE  TO  USING  A  3-TERM  PARTITION  FUNCTION 
FOR  N  AND  A  4-TERM  PARTITION  FUNCTION  FOR  N+  (THE 
IONISATION  ENERGIES  ARE  CORRECTLY  DEPRESSED) 
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FIG. 20  ERRORS  IN  THE  OCCUPATION  NUMBER  OF  THE  iTH  ENERGY 
LEVEL  OF  N  DUE  TO  USING  VACUUM  IONISATION  ENERGIES, 
A  19-TEFIM  PARTITION  FUNCTION  FOR  N  AND  A  4-TERM 
PARTITION  FUNCTION  FOR  N* 
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